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(slight  to  highly  favorable  pressure  gradients).  Results  presented  include 
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EFFECTS  OF  STREAMWISE  PRESSURE  GRADIENTS  ON 


TURBULENCE  STRUCTURE  IN  THE  VISCOUS  LAYER 

D.M.  McEligot  and  H.  Eckelmann 

1.  BACKGROUND 


The  structure  of  a  turbulent  boundary  layer  (TBL)  determines  its  radiated 
noise  and  self  noise  characteristics.  For  most  submersibles ,  if  boundary  layer 
control  is  not  employed  effectively,  the  boundary  layers  near  the  nose  and  near 
hull  arrays  are  turbulent.  Along  a  submersible,  the  streamwise  pressure  gradi¬ 
ent  is  typically  highly  favorable  at  first  (inhibiting  transition  and  stabiliz¬ 
ing  TBLs  after  transition)  and  successively  less  favorable  until  the  location 

of  C  .is  reached;  then  it  is  unfavorable. 
p,mm 

"Strong"  favorable  pressure  gradients  have  been  found  to  reduce  viscous  drag 
by  altering  the  structure  of  the  turbulent  velocity  fluctuations.  Reductions 
in  velocity  fluctuations  and  wall  shear  stress  fluctuations  are  expected  to 
reduce  pressure  fluctuations.  Evidence  for  this  expectation  is  seen  in  models 
of  the  wave- vector-frequency  spectrum  of  wall  pressure  fluctuations  that  vary 
with  the  friction  velocity,  u  =  (t  /p)l/3,  such  as  the  one  by  Chase  [1980]. 

T  W 

With  a  favorable  non-dimensional  pressure  gradient  K  on  the  order  of  -0.03 

P 

an  apparent  "laminarization"  of  a  turbulent  boundary  layer  can  occur  [Narasimha 
and  Sreenivasan,  1979;  McEligot,  1986]. 

Leehey  [1985]  notes  that  the  low  wavenumber  component  of  the  wall  pressure 
spectrum  is  a  potentially  important  source  of  structural  excitation  and  conse¬ 
quent  inward  radiation  of  sound  to  sonar  transducers.  Landahl  [1975J  suggested 
that  unsteady  shear  stresses  should  contribute  significantly  both  to  radiation 
and  to  the  low  wavenumber  component  of  the  wall  pressure  spectrum.  Measurements 
by  Martini,  Leehey  and  Moeller  [1984]  show  that  the  low  wavenumber  component  may 
be  consistent  with  a  shear  stress  source. 

Data  on  the  effects  of  streamwise  mean  pressure  gradients  on  wall  pressure 
fluctuations  are  limited.  For  example,  the  extensive,  critical  review  by 
Willmarth  [1975]  cites  measurements  only  by  Schloemer  [1966],  Bradshaw  [1967], 
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Lim  [1971]  and  Burton  [1973]  for  the  effect  of  streamwise  pressure  gradient  on 

wall  pressure  fluctuations.  The  more  recent  review  of  Simpson,  Ghodbane  and 

McGrath  [1987]  adds  only  airfoil  data  for  adverse  pressure  gradients  by  Hahn 

[1976]  and  Brooks  and  Hodgson  [1981]  to  this  short  list.  Their  own  data  on 

wall  pressure  frequency  spectra  concentrate  on  a  separated  turbulent  boundary 

layer;  they  found  that  t  seemed  to  be  a  better  measure  of  the  stress  to 

max 

use  to  scale  p'  than  t  (however,  for  a  favorable  pressure  gradient, 

w 

T  -  T  )  . 

max  w 

Blake  [1986]  compared  spectra  from  measurements  in  adverse  pressure  gradients 
without  successfully  correlating  them.  He  noted  that  the  spectrum  shape  appar¬ 
ently  depends  on  the  upstream  history  of  the  boundary  layer  development,  since 
there  is  no  correlation  of  the  results  with  Reynolds  number,  wall  shear  stress 
or  pressure  gradient. 

There  are  three  well  known  models  of  wall  pressure  spectra  describing  the 
turbulent  boundary  layer  flow  noise;  they  are  by  Corcos  [1963,  1967],  Chase 
[1980]  and  Ff owes- Williams  [1982].  An  excellent  review  and  assessment  of  the 
engineering  applicability  of  the  above  three  models  was  given  by  Hwang  and  Geib 
[1983].  They  examined  the  wavevector-f requency  spectrum  of  the  turbulent 
boundary  layer  wall  pressure  in  the  intermediate  and  high  frequency  range. 

Their  intercomparisons  showed  that  both  the  Corcos  and  Chase  models  are  suitable 
for  engineering  applications;  however,  computation  with  the  Ffowcs-Williams 
model  could  not  be  pursued  due  to  the  still-unknown  functions  required  in  that 
model.  In  general,  the  Chase  model  gave  better  results  and  the  spectrum 
corresponding  to  q^  is  sufficiently  accurate  for  engineering  applications. 
Engineering  application  of  this  modified  spectrum  to  establish  the  flow  noise 
over  a  submersible  may  be  a  feasible  approach  coupled  with  other  correction 
factors  due  to  pressure  gradients,  surface  curvature,  etc. 

For  negligible  pressure  gradients,  theoretical  models  due  to  Kraichman 
[1956],  Meng  and  Lovett  [1985]  and  others  have  demonstrated  relationships 
between  the  wall-pressure  frequency  spectrum  and  the  integral  of  the  mean 
square  velocity  in  the  boundary  layer.  From  the  scaling  analysis  of  Meng  and 
Lovett,  one  finds  that  for  incompressible  flow  the  wall  pressure  spectrum  may 
be  expressed  as 
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p(0,k,w)  =  pk  J  e  [JJJ  u  e  d  xdtjdy 

2  (  2ir) 3  0  -oo 

where  tc  is  the  magnitude  of  the  two-dimensional  vector  wavenumber  in  the 
plane  of  the  wall  and  <*>  is  the  frequency.  This  relation  is  the  basis  of  a 
proposal  by  Meng  to  determine  the  wall  pressure  spectrum  via  LDA  measurements. 
His  order-of-magnitude  analysis  indicates  that  high  wavenumber  variations  are 
due  to  fluctuations  in  the  linear  layer  (y+<5),  convective  peak  variations 
from  the  viscous  layer  (y+<30)  and  low  wavenumber  contributions  from  the  outer 
region.  However,  with  favorable  pressure  gradients  these  regions  are  thickened 
so  their  related  wavenumber  spectral  regions  are  shifted,  e.g.,  the  viscous 
layer  can  begin  to  affect  the  low  wavenumber  component. 

For  significant  transverse  curvature,  as  in  a  small  diameter  towed-array, 
Chase  [1981]  examined  the  possible  scaling  of  the  axisymmetric  component  of 
spectral  density  of  wall  pressure  with  a  low- wavenumber  expansion  in  terms  of 
fluctuating  Reynolds  stresses  as  sources.  A  model  of  the  sources  was  formed 
based  on  similarity  near  the  wall  (e.g.,  fixed  mixing  length  function).  For 
low  wavenumbers  he  developed  a  relationship  of  the  form 


P°  (k,u)  =  C  PV  fn(k,w} 

o  2  T 

where  uT  is  the  friction  velocity  /t w/p . 

Panton  and  Linebarger  [1974]  derived  wall-pressure  wavenumber  spectra  for 
zero-pressure  gradient  and  adverse-pressure-gradient  equilibrium  boundary  layers 
based  on  Coles'  laws  of  the  wall  and  wake  and  on  turbulence  structure  being 
invariant  with  pressure  gradient.  There  appear  to  be  no  models  of  the  wall- 
pressure  wavenumber  spectra  that  take  into  account  the  effects  of  favorable 
streamwise  pressure  gradients  on  the  turbulence  contributions  to  the  source 
terms . 

Handler,  Hansen,  Sakell,  Orzag  and  Bullister  [1984]  obtained  predictions 
of  the  wall  pressure  field  p(x,z,t)  for  fully-developed  flow  in  a  channel  by 
direct  solution  of  the  time-dependent  Navier-Stokes  equations.  Presumeably, 
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the  comparable  calculations  by  Kim,  Moin  and  Moser  [1987]  could  be  analyzed  to 
provide  the  same  sort  of  information;  for  tu^/w  =  10,  this  case  required 
0(2  CPU  weeks)  with  a  CRAY-XMP  computer.  How  to  extend  these  techniques  to 
turbulent  boundary  layers  with  significant  pressure  gradients  is  not  a  trivial 
question  (but  the  work  herein  should  provide  useful  tests  of  some  aspects  of 
such  predictions). 

Hydrodynamic  boundary  layer  predictions  can  be  used  to  predict  the  local 
boundary  layer  characteristics  with  pressure  gradients,  streamwise  curvature, 
etc.,  so,  conceptually,  models  such  as  the  ones  by  Chase  [1980,  1981]  and 
Panton  and  Linebarger  [1974]  could  be  used  to  predict  wall  pressure  spectra. 

The  question  is  whether  the  turbulence  structure  is  invariant  under  these 
phenomena  as  implied  in  these  models.  Measurements  of  the  fluctuating  wall 
shear  stress  in  favorable  pressure  gradients  by  Chambers,  Murphy  and  McEligot 
[1983]  indicate  it  probably  is  not;  for  example,  the  non-dimensional  bursting 
frequency  varies  with  pressure  gradient. 

Only  a  few  experiments  have  measured  the  instantaneous  turbulence  structure 
simultaneously  with  the  wall  pressure  signal,  and  these  have  been  in  zero 
pressure  gradients  or  fully  developed  flow.  Zakkay,  Barra  and  Hozumi  [1979], 
Thomas  and  Bull  [1983]  and  Kobashi,  Komoda  and  Ichijo  [1984]  essentially  used 
wind  tunnels  and  Dinkelacker  and  Langeheineken  [1982;  bangeheineken ,  1981]  at 
the  Max  Planck  Institut  fur  Stromungsf orschung  conducted  their  measurements  in 
fully-developed  air  flow  in  a  circular  tube.  All  confirmed  that  the  wall 
pressure  fluctuations  were  related  to  the  turbulent  velocity  fluctuations. 
(Since  different  detection  criteria  were  used  these  investigations  cannot  be 
compared  directly.) 

Thomas  and  Bull  [1983]  measured  instantaneovs  wall  pressure  and  wall  shear 
stress  simultaneously,  then  by  conditionally  sampling  and  averaging  they  were 
able  to  examine  and  compare  their  characteristic  signatures.  Thomas  and  Bull 
and  the  other  investigators  showed  that  the  wall  pressure  signal  was  related  to 
the  bursting  phenomena  [Blackwelder  and  Kaplan,  1976]  of  turbulent  flow. 
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2.  PREVIOUS  WORK 


As  noted  above,  it  is  possible  to  predict  the  fluctuating  pressure  field 
if  the  turbulence  structure  of  the  fluctuating  velocity  field  is  known.  This 
section  summarizes  the  studies  of  the  turbulence  structure,  aimed  towards  con¬ 
centration  on  flows  with  favorable  streamwise  pressure  gradients. 

In  an  internal  flow,  a  favorable  streamwise  pressure  gradient  is  induced 
by  a  convergence  -  either  lateral  or  transverse  -  leading  to  a  spatial  accelera 
tion  of  the  flow.  For  external  flows,  the  body  shape  near  the  nose  or  leading 
edge  induces  a  favorable  pressure  gradient  and,  via  Euler's  equation,  accelera¬ 
tion  of  the  freestream  flow. 

While  the  details  of  many  accelerating  flows  have  been  extensively  studied 
[Narasimha  and  Sreenivasan,  1979;  Murphy,  1979),  only  McEligot  and  colleagues 
appear  to  have  utilized  lateral  convergence  between  parallel  plates  to  examine 
turbulent  flows  with  positive  streamwise  pressure  gradients  [Murphy,  Chambers 
and  McEligot,  1983;  Chambers,  Murphy  and  McEligot,  1983).  In  the  idealized 
problem  (Fig.  2.1)  the  velocity  is  uniform  at  an  entry  at  a  large  distance. 

The  flow  develops  and  may  approach  a  fully  established  profile.  Near  the  exit 
of  the  convergence,  acceleration  dominates  the  flow  in  the  center  of  the  duct, 
causing  a  flattening  of  the  profile.  The  profile  takes  the  appearance  of  an 
external  boundary  layer  flow,  with  a  large  central  core  flanked  on  both  sides 
with  thin  boundary  layers.  The  boundary  layers  become  successively  thinner  as 
acceleration  continues.  Depending  on  the  entrance  length  and  flow  rate,  the 
regime  may  be  laminar  or  turbulent,  or  an  initially  turbulent  flow  may  approach 
a  laminar  flow  owing  to  the  stabilizing  effects  of  acceleration,  i.e., 

' laminarization '  may  occur. 

Several  authors  have  proposed  various  non-dimensional  parameters  involving 
either  velocity,  pressure  or  shear-stress  gradients  which  purportedly  serve  as 
guides  to  the  likelihood  of  laminarization.  Kline  et  al.  [1967]  noticed  that 
turbulent  bursts  ceased  when  an  acceleration  parameter  K  =  (v/u2)du  /dx 

\j  oo  oo 

reached  a  critical  value  of  3.7  x  10  .  Jones  and  Launder  [1972]  found  that 

the  peak  Reynolds  stress  due  to  turbulence  became  a  progressively  smaller 

fraction  of  the  wall  shear  stress  as  K  increased.  They  indicated  that  a  turbu 

u 
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lent  flow  could  no  longer  exist  when  K  exceeded  a  critical  value  somewhere 

u 

between  2.5  x  10  and  3.5  x  10  .  Blackwelder  and  Kovasznay  11972}  found 

that,  when  exceeded  the  above  critical  values,  the  viscous-sublayer  thick¬ 
ness  increased,  the  skin-friction  coefficient  at  the  wall  decreased,  and 
large  departures  from  the  turbulent  logarithmic  law  of  the  wall  occurred. 

Patel  [1965]  also  noted  major  departures  from  the  logarithmic  law,  but 

correlated  his  results  in  terms  of  a  pressure-gradient  parameter 

Kp  =  (v/pu^)dp/dx. Patel  and  Head  [1968]  suggested  0.018  for  a  critical 

value  of  K  .  They  also  suggested  that  a  more  appropriate  laminarization 

p 

parameter  might  be  defined  by  replacing  the  pressure  gradient  with  the  shear 

stress  gradient  at  the  wall  Ox/Syl^;  we  shall  call  it  K  .  Since  the 

gradient  3-r/3y  is  evaluated  at  the  wall,  where  the  convective  terms  in 

the  boundary- layer  momentum  equation  are  zero,  dp/dx  =  (d-r/dy)  for  flow 

w 

between  parallel  plates,  whether  fully  developed  or  accelerating;  for  fully 

developed  pipe  flow  dp/dx  =  2Mt(y)/My.  For  transitional  pipe  flow  Patel 

and  Head  found  that  with  the  present  nomenclature  K  =  2K  =  0.018,  whereas 

P  t 

for  transitional  flow  between  parallel  plates  they  found  that  =  K  =  0.009. 

On  the  basis  of  consistency,  then,  they  concluded  that  was  a  better  guide 

to  the  possibility  of  laminarization  and  that  its  critical  value  was  0.009 

or  so.  They  also  suggested  that  the  downstream  distance  required  to  reverse 

transition  completely  must  be  dependent  upon  the  Reynolds  number  based  on 

momentum  thickness,  R  =  u  6  /v.  Narasimha  and  Sreenivasan  [1973]  pointed 

2  co  2 

out  that  the  critical  value  of  should  have  been  0.004,  owing  to  algebraic 

errors.  They  concluded  that  laminarization  was  due  to  domination  of  pressure 

gradients  over  the  Reynolds  stresses  and,  therefore,  the  parameter  K  incorpo- 

P 

rating  the  pressure  gradient  was  recommended.  In  addition,  they  pointed 
out  that  parameters  such  as  K  ,  K  and  K  cannot  be  universally  useful  for 

up  T 

all  levels  of  acceleration.  When  the  acceleration  is  very  low,  as  in  a 
nearly  fully  developed  duct  flow,  the  flow  and  its  state  of  turbulence 
revert  to  their  equilibrium  values  and  the  local  Reynolds  number  then  pro¬ 
vides  the  best  guide  to  the  possibility  of  laminarization. 

From  measurement  of  the  local  Nusselt  number  in  strongly  heated  gas  flow 
in  a  circular  tube,  McEligot  [1963]  detected  premature  transition  from  turbu¬ 
lent  flow  towards  laminar  flow.  He  proposed  the  explanation  to  be  a  thickening 
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the  viscous  sublayer  and  correlated  the  flow  regimes  in  terms  of  a  non- 
dimensional  heating  rate  and  inlet  Reynolds  number.  The  correlation  was 
later  shown  by  McEligot,  Coon  and  Perkins  [1970]  to  be  related  to  an  acceler¬ 
ation  parameter  Kv  =  (v/V2)dV/dx  based  on  the  bulk  velocity.  This  para¬ 
meter  K  is  a  version  of  K  appropriate  for  internal  flows.  Another 
v  u 

Reynolds- number  criterion  was  proposed  by  Bradshaw  [1969],  who  postulated 

that  a  turbulent  flow  would  become  directly  dependent  upon  viscosity,  i.e., 

laminarizing,  when  the  energy- containing  and  the  dissipation  ranges  of  eddy 

size  overlap.  This  led  to  a  critical  ’eddy’  Reynolds  number  R  ..  =  u  L/v  =  12, 

eddy  t 

where  L  is  the  dissipation-length  parameter  u^e  and  e  is  the  dissipation  rate. 

Badri  Narayanan  and  Ramjee  [1969]  attempted  to  meld  these  diverse  observa¬ 
tions  into  a  unified  structure  by  describing  reverse  transition  in  three  con¬ 
secutive  states:  (i)  the  disappearance  of  the  turbulent  bursts  near  the  wall  at 

K  ^  3  x  10  *;  (ii)  breakdown  of  the  law  of  the  wall  at  K  =  0.02;  and  (iii) 
u  p 

the  decay  of  turbulence  intensity  starting  at  R  =  300-400.  While  at  first 

2 

glance  this  description  seems  to  accommodate  all  the  observations  reported 

above,  stage  (iii)  strongly  implies  that  the  completion  of  laminarization  is 

controlled  by  a  critical  value  of  R  ,  whereas  Launder  and  Jones  [1969] 

2 

suggested  that  was  the  critical  parameter.  In  summary,  there  still  remain 
significant  uncertainties  not  only  in  the  precise  values  of  the  dimensionless 
parameters  to  employ  as  criteria,  but  also  as  to  which  parameter  is  the  best  to 
use  in  predicting  the  occurrence  of  laminarization. 

For  single  phase  boundary  layers  on  flat  plates  with  zero  pressure  gradi¬ 
ent  or  in  fully-developed  duct  flow,  a  large  number  of  studies  have  detailed 
the  turbulence  structure  by  flow  visualization  and  hot-wire  or  hot-film 
anemometry  [Corino  and  Brodkey,  1969;  Kim,  Kline  and  Reynolds,  1971;  Eckelmann, 
1974;  Johansson  and  Alfredsson,  1982;  etc.].  At  the  Max  Planck  Institut  fur 
Stromungsforschung,  Eckelmann  [1974]  and  his  colleagues  [Wallace,  Brodkey  and 
Eckelmann,  1977;  Randolph,  1983;  Haritonidis  and  Johansson,  1985,  etc.]  have 
developed  techniques  for  structure  measurements  in  their  unique  oil  channel, 
designed  originally  by  Reichardt  to  give  excellent  spatial  resolution.  The  use 
of  oil  also  leads  to  higher  temporal  resolution  of  fluctuating  components  than 
in  the  common  fluids,  air  and  water. 
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The  bursting  phenomenon  in  the  viscous  layer  has  been  described  well  by 
Corino  and  Brodkey  11969],  Kline  et  al.  {1967],  Often  and  Kline  [1974]  and 
others,  as  well  as  in  reviews  by  Willmarth  [1975]  and  Cantwell  [1981],  so  only 
the  literature  pertinent  to  the  present  study  will  be  reviewed  here.  With  Kim, 
Kline  and  Reynolds  [1971]  and  Blackwelder  and  Kaplan  [1976],  we  adopt  the  term 
•bursting'  or  ’bursting  phenomenon*  for  the  overall  process.  For  convenience 
in  description,  the  phases  of  the  process  will  be  termed  (a)  a  deceleration  or 
ejection,  (b)  a  rapid  acceleration  or  sweep,  and  (c)  a  relatively  quiescent 
process,  or  more  gradual  deceleration  following  the  sweep.  Blackwelder  and 
Kaplan  [1976]  have  demonstrated  one  version  of  the  typical  bursting  signal  by 
conditional  averaging  of  the  fluctuation  in  the  streamwise  velocity  component. 
(Hereinafter  we  will  usually  refer  to  their  paper  as  BK.)  With  calculations 
based  on  a  theoretical  model  for  coherent  structures  in  wall  turbulence,  Landahl 
[1980]  has  found  qualitative  agreement  with  their  conditionally  sampled  data. 

From  examination  of  the  simultaneous  measurements  of  u(t)  and  3u(t)/3y, 
Eckelmann  [1974]  showed  that  the  fluctuating  shear-stress  signal  from  a  wall 
shear- stress  sensor  correlates  with  velocity  fluctuations  measured  with  hot 
wires  positioned  in  the  flow  at  dimensionless  distances  from  the  wall  of  up  to 
y+  =  25.  Eckelmann' s  correlation  indicates  that  the  temporal  history  of 
either  the  shear-stress  signal  or  the  fluctuating  streamwise  velocity  component 
(u)  should  provide  a  reasonable  picture  of  turbulent  events  in  the  viscous 
layer,  0  <y+  <30  in  unaccelerated  flows.  In  addition,  Brown  and  Thomas 
[1977]  have  shown  that  the  low-frequency  fluctuations  of  a  wall-shear  sensor 
can  be  correlated  with  the  low-frequency  fluctuations  of  a  velocity  sensor  for 
the  range  0.05  <y/6  <0.75  in  a  turbulent  boundary  layer.  Eckelmann  [1980] 
notes  that  the  detailed  similarity  between  the  two  signals  decreases  with 
distance  from  the  wall,  so  one  must  investigate  whether  one  sees  the  same  burst 
at  y+  =  15,  the  triggering  location  of  BK,  and  at  the  wall. 


A  second  question  concerns  the  appropriate  scaling  of  the  bursting 

frequency  f  in  non-accelerated  flows.  As  discussed  below,  some  investigators 

have  suggested  non-dimensionalization  with  respect  to  a  characteristic  time  for 

the  wall,  or  inner  region,  v/Cu^)2;  others  advocate  a  time  characteristic  of  the 

outer  region,  6/U  .  For  a  fully  developed  internal  flow,  the  two  resulting 
00 
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The  question  has  been  one  of  whether  either  scaling  will  give  a  non-dimen¬ 
sional  time  or  frequency  which  is  independent  of  Reynolds  number.  The  question 
of  scaling  is  an  important  one  for  a  better  understanding  of  turbulent  behavior 
because  inner  scaling  implies  that  the  lift-up  of  the  low-speed  streak  and  the 
subsequent  ejection  [Offen  and  Kline,  1974]  are  derived  from  a  wall-flow  insta¬ 
bility;  theoretical  models  showing  such  instabilities  have  been  proposed  [Black 
1968],  (However,  these  instabilities  cannot  be  unambiguously  accepted  as  the 
causes  of  bursting  because  the  instability  models  usually  require  [Black  1968] 
that  the  streamwise  velocity  of  the  ejected  fluid  be  greater  than  the  local 
mean  velocity  corresponding  to  the  y+  position  of  the  ejection,  rather  than 
the  velocity  deficit  actually  observed.)  Supporting  evidence  for  inner  scaling 
is  offered  by  the  flow-visualization  data  of  Offen  and  Kline  [1974],  who  found 
that  sweeps,  low- speed-streak  lift-ups,  and  bursts  were  intimately  connected; 
indeed,  ’sweeps  initiate  bursts',  and  'sweeps  originate  in  the  inner  layer’. 


Outer  scaling  implies  that  the  bursting  process  is  imposed  by  large  dis¬ 
turbances  already  present  in  the  outer  flow.  However,  even  Rao,  Narasimha  and 
Badri  Narayanan  [1971],  early  proponents  of  outer  scaling,  admitted  that  the 
’origin  of  the  bursting  phenomenon  cannot  be  traced  directly ...  to  the  outer 
interface’.  Outer  scaling  has  its  theoretical  support,  too.  The  analysis  by 
Schubert  and  Corcos  [1967]  attempted  to  view  the  inner  layer  as  driven  by  the 
outer  flow,  but  failed  to  predict  Reynolds  stresses  of  the  right  order.  Brown 
and  Thomas  [1977]  inferred  from  their  measurements  that  the  low-frequency 
fluctuations  of  the  wall  shear  stress  are  possibly  driven  by  large-scale 


structures . 


Attempts  to  resolve  the  scaling  issues  experimentally  have  yielded  mixed 

results.  Black  [1966]  found  that  the  bursting  frequencies  scaled  with  wall 

2 

variables  so  that  (u^)  /vf^  ^  110.  Kline  et  al.  [1967]  reported  essentially 

the  same  finding  for  a  boundary- layer  flow,  as  did  Corino  and  Brodkey  [1969]  for 
fully  developed  flow  in  a  pipe.  Sharma  and  Willmarth  [1980]  found  the  non- 
dimensional  bursting  frequency  to  be  approximately  constant  over  the  range 
1  <  y+  <  17  when  based  on  wall  variables  (this  observation  provides  further 
support  for  the  use  of  u  wall-mounted  sensor  to  investigate  the  bursting  process). 
In  contrast,  Rao  et  al.  [1971]  suggested  that  outer  scaling  resulted  in 
Reynolds  number  independence  for  unaccelerated  flows,  and  that  U  /6fL  =  5 

oo  b 

or  U  /&*f,  =  32,  where  S*  is  the  displacement  thickness. 

oo  t) 

A  variety  of  techniques  has  evolved  to  identify  and  measure  turbulence  and 
its  coherent  structures.  Among  the  most  popular  are  the  VITA  technique 
[Blackwelder  and  Kaplan,  1976]  and  quadrant  splitting  [Wallace,  Eckelmann  and 
Brodkey,  1972].  There  has  been  concern  about  the  various  detection  techniques 
and  their  consequent  conditional  averages  [Johansson  and  Alfredsson,  1982]  and 
several  studies  with  negligible  pressure  gradients  have  been  devoted  to  their 
comparison  [Subramanian  et  at.,  1982;  Alfredsson  and  Johansson,  1984;  Luchik 
and  Tiederman,  1987].  Luchik  and  Tiederman  found  all  they  examined  to  yield 
the  same  average  time  between  bursts  but  different  techniques  showed  different 
advantages . 

Additionally,  the  finding  of  Rao  et  al.,  was  based  upon  wind-tunnel 
measurements  in  which  a  single  hot-wire  signal  was  differentiated  and  a  dis¬ 
crimination  level  varied  to  obtain  a  maximum  f.  for  each  run.  (It  was  stated 

D 

that,  as  the  discrimination  level  was  set  either  too  low  or  too  high,  f. 

b 

would  decrease.)  Offen  and  Kline  [1974]  indicate  that  such  burst  detection  by 
a  single  sensor  using  autocorrelations  has  little  chance  of  success,  except  very 
near  the  wall,  because  the  individual  bursts  vary  greatly  in  size.  Thus 
detection  techniques  based  upon  spatial  rather  than  temporal  coherence  are 
preferred.  Flow-visualization  studies,  such  as  those  by  Kline  et  al.  [1967] 
and  Corino  and  Brodkey  [1969],  emphasize  the  more  easily  detected  spatial 
coherence  of  bursts.  A  fixed  point  sensor  can  only  measure  the  less  organized 
temporal  behavior  of  the  flow.  Using  the  same  apparatus  as  Kline  et  al.,  Kim 
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et  al.  [1971]  determined  bursting  frequencies  for  two  flows  from  visualization 
studies  and  from  the  disputed  autocorrelation  method  with  a  hot  wire  at 
y+  =  14.  Because  the  Reynolds  numbers  of  the  two  flows  were  nearly  equal, 
the  results  could  be  interpreted  as  supporting  either  inner  or  outer  scaling. 

Strickland  and  Simpson  [1975]  measured  bursting  frequencies  in  two  ways: 

(i)  the  time  between  peaks  in  the  autocorrelation  of  a  wall-mounted  shear- 

stress  sensor  was  taken  as  a  measure  of  the  time  between  bursts,  giving  f  , 

and  (ii)  f.  was  taken  to  be  the  frequency  at  the  maximum  in  the  first 
bs 

moment  of  the  curve  of  power-spectral-density  versus  frequency.  They  found 
agreement  of  the  two  frequencies  to  within  15%.  While  U  /if.  ranged 

CD  D 

from  7  to  14,  in  fair  agreement  with  Rao  et  al.,  most  of  the  flows  examined 
were  decelerating  flows.  For  the  two  runs  where  the  flow  was  mildly 
accelerated,  U  /if.  was  approximately  11,  twice  the  value  recommended 

oo  D 

by  Rao  et  al.  The  Reynolds  numbers  of  the  two  flows  were  close  enough 
(Re  =2240  and  3520)  that,  had  a  Reynolds-number  dependence  existed,  it 
would  not  have  been  easily  apparent.  Indeed,  in  related  work  Strickland 
[1973]  reported  a  weak  Reynolds  number  dependence. 

Wallace,  Brodkey  and  Eckelmann  [1977]  attempted  to  recognize  bursts  using 
a  pattern-recognition  technique  based  upon  the  time  derivatives  of  hot-wire 
velocity  sensors.  They  found  that  (Wif^  ranged  from  5  at  y+  =  3  to 
2  at  y+  =  300.  Although  the  agreement  with  Rao  et  al.  was  apparently  good, 
only  a  single  flow  was  investigated,  so  that  Reynolds-number  independence 
was  not  substantiated.  Blackwelder  and  Kaplan  [1976],  using  the  variable- 
interval  time-averaging  (VITA)  technique,  also  reported  that  IWif^  *  2.8, 
but,  again,  essentially  only  one  Reynolds  number  was  examined. 

With  a  direct  numerical  simulation,  Kim  and  Spalart  [1987]  examined  the 
scaling  of  the  bursting  frequency  for  momentum  Reynolds  numbers,  300  <  Re-.  < 

1410.  They  concluded  that  the  frequency  scaled  by  inner  variables  is  essentially 
independent  of  Reynolds  number. 

Quantitative  turbulence  structure  measurements  with  pressure  gradients  are 
sparse.  In  general,  they  are  limited  to  data  with  a  single  sensor  (no  statis- 
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tics)  due  to  probe  size  concerns  using  air  or  water  as  the  fluid.  Jones  and 
Launder  [1972]  found  that  the  peak  Reynolds  stress  due  to  turbulence  became  a 
progressively  smaller  fraction  of  the  wall  shear  stress  as  K  =  (v/u*)du  /dx 

VJ  oo  oo 

increased.  In  a  study  of  large  scale  motion  Blackwelder  and  Kovasznay  [1972] 
found  that,  when  Ku  exceeded  "critical"  values,  the  viscous-sublayer  thickness 
increased,  the  skin-friction  coefficient  at  the  wall  decreased,  and  large 
departures  from  the  turbulent  logarithmic  law  of  the  wall  occurred;  to  some 
extent  these  aspects  can  be  predicted  using  advanced  (and  some  simple)  turbu¬ 
lence  models  [Rodi,  1980;  McEligot,  1986].  Simpson  [1979]  observed  that  the 
spacing  of  streamwise  vortices  in  the  wall  region  seems  to  increase  with  K. 
Finnicum  and  Hanratty  [1987]  measured  wall  shear  stresses  in  pulsating  flows 
and  compared  them  to  steady  flow.  Recently  Tiederman  and  colleagues  have  begun 
the  measurement  of  turbulence  structure  with  adverse  pressure  gradients 
[White  and  Tiederman,  1990]. 

Analyses  predicting  the  effects  of  streamwise  pressure  gradient  on  turbu¬ 
lent  boundary  layers  are  now  available  from  the  work  of  Spalart  [1986]  and  of 
Finnicum  and  Hanratty  [1988].  Spalart  conducted  direct  numerical  simulations 
of  sink-flow  boundary  layers  with  acceleration  parameters  between  1.5  x  10 
and  3.0  x  10  6 .  He  solved  the  three-dimensional,  time-dependent  Navier-Stokes 
equations  using  a  spectral  method.  Predicted  effects  of  the  favorable  pressure 
gradients  were  to  extend  and  displace  the  logarithmic  layer  and  to  alter  the 
energy  balance  of  turbulence  near  the  edge  of  the  boundary  layers.  Relaminar- 
ization  was  predicted  for  =  3.0  x  10  6.  Reasonable  agreement  was  found  with 
the  spatial  spectra  deduced  by  Jones  and  Launder  [1972]  from  their  time  spectra 
data  at  =  1.5  x  10  6  except  that  Spalart' s  computed  spectra  did  not  collapse 
versus  y.  In  general  his  results  suggested  that  the  logarithmic  behavior  of  the 
mean  velocity  is  more  universal  than  the  linear  behavior  of  the  mixing  length. 

Finnicum  and  Hanratty  applied  the  2-1/2D  computational  model  of  Nikolaides 
[1984]  with  bursting  periods  and  streak  spacing  taken  to  vary  as  the  ratio  be¬ 
tween  the  wall  shear  stress  and  the  mean  shear  stress  at  the  assumed  edge  of 

the  viscous  layer  (y+) .  Empirical  results  were  used  for  skin  friction  coef- 

o 

ficients.  Phase  relations  at  y^  were  assumed  to  be  unchanged  by  the  pressure 
gradient.  For  a  range  K  =0,  2.0  x  10  6  ("moderate")  and  2.8  x  10  6  ("close 
to  relaminarization")  they  predicted  mean  statistics  of  u  and  v  and  appropriate 
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terms  for  kinetic  energy  balances.  Good  agreement  was  seen  with  the  direct 
numerical  simulations  of  Spalart  [1986). 

3.  STUDIES  BY  PRINCIPAL  INVESTIGATOR  AND  COLLEAGUES 


For  laterally  converging  laminar  flow  between  parallel  plates,  Murphy, 

Coxon  and  McEligot  [1978]  extended  the  numerical  approach  of  Bankston  and 
McEligot  [1970],  as  well  as  developing  a  similarity  solution.  Murphy,  Chambers 
and  McEligot  [1983]  then  applied  the  numerical  program  to  the  conditions  of  an 
experiment  they  conducted  with  turbulent  and  apparently  laminarizing  air  flow 
converging  between  two  plates  separated  by  1.3  cm  (1/2  in.).  The  test  section 
had  a  sharp-edged  entrance  to  induce  a  transition  to  turbulent  flow  and  side 
rails,  whose  angle  could  be  adjusted  to  vary  the  acceleration  parameter  and 
Reynolds  number  independently  (Figure  3.1).  With  streamwise  acceleration  the 
flow  was  continuously  developing  through  the  entry  region  and  downstream. 
Measurements  were  restricted  to  streamwise  pressure  distributions  and  wall 
shear  stress,  since  the  spacing  necessary  was  too  small  to  introduce  probes 
into  the  flow  to  determine  the  velocity  distribution.  By  comparison  between 
data  and  predictions,  they  found  substantial  effects  to  occur  with  convergence 
angles  as  small  as  4°  and  showed  that  the  results  could  be  explained  by  a 
thickening  of  the  viscous  layer. 

In  the  same  apparatus  Chambers,  Murphy  and  McEligot  [1983]  applied  condi¬ 
tional  sampling,  by  the  VITA  technique  of  Blackwelder  and  Kaplan  [1976],  to  the 
instantaneous  measurements  of  wall  shear  stress  (Figure  3.2).  They  found  the 
typical  burst  pattern,  or  conditionally-averaged  time  history  of  the  wall  shear 
stress,  resembled  the  time  history  of  the  streamwise  velocity  component  measured 
by  Blackwelder  and  Kaplan  at  y+  =  15.  Families  of  conditionally-averaged 
time  histories  were  selected  at  approximately  equal  Reynolds  numbers  and  varying 
acceleration  parameters  to  compare  with  the  results  for  fully  developed  flows. 
Figure  3.2  presents  typical  results  at  a  Reynolds  number  of  about  12,000.  Addi¬ 
tional  families  were  compared  at  Re  =  8,000,  18,000  and  23,000  with  the  same 
qualitative  conclusions. 

At  a  given  Reynolds  number,  there  is  no  large  variation  in  the  shape  of  the 
time- averaged  signals  until  the  acceleration  parameter  reaches  L  ?  2  x  10  *. 
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That  is,  the  sweep  pattern  is  essentially  invariant  until  Ky  approaches  this 
value  even  though  the  non-dimensional  bursting  frequency  begins  to  drop  con¬ 
siderably  for  Ky  >  1  x  10  6.  At  higher  values  than  Ky  -  2  x  10  6  the  recovery 
or  decay  of  the  sweep  phase  becomes  slower  than  in  the  fully  developed  or  moder¬ 
ately  converging  runs. 

For  fully  developed  flows  (convergence  angle  =  0)  wall  scaling  of  the 
bursting  frequency  was  found  to  be  essentially  independent  of  Reynolds  number. 
The  effects  of  laterally  converging  flows  upon  bursting  frequencies  are  apparent 
in  Figure  3.3.  No  matter  how  they  are  scaled,  all  the  dimensionless  frequen¬ 
cies  suggest  a  significant  decrease  in  f.  or  increase  in  its  inverse,  the 

D 

average  period,  as  Ky  increases. 

Focusing  upon  the  inner  scaling  of  ffa,  shown  to  be  the  most  appropriate 

scaling  in  the  case  of  fully  developed  flows,  one  finds  that  as  Ky  decreases 

towards  10  7  the  bursting  frequency  approaches  the  values  observed  for 

—  6  2 

fully  developed  turbulent  flows.  At  Ky  =  1  x  10  ,  u^/vf^is  300,  nearly 
70%  higher  than  the  value  at  Ky  =  0;  at  Ky=  2  x  10  6  one  finds  a  threefold 
increase,  and  Ky  -  3.4  x  10  6  indicates  a  fifteenfold  increase,  with  the  curve 
beginning  to  tend  asymptotically  towards  infinity.  This  tendency  implies  that 
the  flow  would  be  effectively  laminar  for  Ky  larger  than  about  3.4  x  10  *. 

This  finding  is  in  accord  with  the  data  of  Murphy,  Chambers  and  McEligot  [1983] 
in  the  same  test  section  and  with  the  data  of  Moretti  and  Kays  [1965],  Kline  et 
al.  [1967],  Blackwelder  and  Kovasznay  [1972]  and  Jones  and  Launder  [1972]  in 
different  geometries  using  other  techniques.  For  laterally  converging  flows, 
the  measurements  of  bursting  characteristics  essentially  confirmed  the  deduc¬ 
tions  of  Murphy,  Chambers  and  McEligot  [1983];  the  apparent  modification  of 
turbulent  behavior  inferred  from  comparison  of  predictions  and  pressure 
measurements  corresponded  with  a  reduction  in  the  non-dimensional  bursting 
frequency. 


For  this  flow  problem,  the  results  of  Murphy,  Chambers  and  McEligot  can  be 
summarized  as  follows:  1)  for  Ky  >  4  x  10  6,  laminar  predictions  are 
adequate  for  pressure  distributions  and  friction  factors,  2)  for  Ky  <  10  7, 
the  simple  van  Driest  [1956]  model  provides  reasonable  predictions,  3)  the  temporal 
wall  shear  signal  was  consistent  with  these  observations,  and  4)  as  Ky 
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increased,  the  non-dimensional  bursting  frequency  decreased  but  the  condi¬ 
tionally-averaged  signal  remained  approximately  constant  until  large  K^. 

In  summary  there  is  evidence  that  the  turbulence  structure  is  modified  by 
streamwise  pressure  gradients  and  there  is  a  need  to  determine  which  features 
remain  relatively  invariant  and  to  quantify  the  variation  of  the  others  in 
order  to  provide  the  turbulence  structure  information  needed  in  hydroacoustic 
models  accounting  for  their  effects. 


4.  GOAL  AND  OBJECTIVES 


The  ultimate  Roal  of  this  study  is  to  develop  practical  models  of  wall 
pressure  spectra  in  turbulent  boundary  layers  where  significant  drag  reducing 
phenomena  may  alter  the  turbulence  structure  that  is  the  basis  of  such  models. 

The  next  steps  are  (1)  to  analyze  the  turbulence  structure  for  measurements  with 
favorable  streamwise  pressure  gradients,  (2)  to  compare  the  results  to  data  with 
drag- reducing  additives,  and  (3)  to  derive  a  model  for  the  effects  of  these  phe¬ 
nomena  on  wall  pressure  spectra  based  on  the  significant  variations  in  turbulence 
structure  deduced  in  steps  (1)  and  (2).  The  present  objective  is  to  accomplish 
the  first  step. 


5.  APPARATUS 


The  experiments  performed  utilized  the  familiar  oil  channel  at  the  Max  Planck 
Institut  fur  Stromungsf orschung  (Figure  5.1a).  This  unique  facility  was  designed 
originally  by  Reichardt  to  give  high  spatial  and  temporal  resolution  in  the 
viscous  layer  and  has  since  been  developed  by  Eckelmann  [1974]  and  his  colleagues 
[Wallace,  Eckelmann  and  Brodkey,  1972;  Wallace,  Brodkey  and  Eckelmann,  1977; 
Blackwelder  and  Eckelmann,  1978;  Randolph,  1983;  etc.] 

The  oil  channel  and  procedures  employed  in  its  use  have  been  described  in 
detail  by  Eckelmann  [1970,  1974].  In  this  section  we  will  concentrate  on  the 
differences  in  apparatus  and  techniques  applied  for  the  present  study;  only 
a  brief  overview  will  be  given  for  the  aspects  which  are  common  with  the 
earlier  work. 
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The  channel  is  22  cm  wide  and  1  m  deep  with  an  8  m  length.  In  contrast  to 
the  1.3  cm  spacing  in  the  experiment  of  Murphy,  Chambers  and  McEligot  [1983], 
the  plates  forming  this  oil  channel  have  much  larger  separation,  so  area  block¬ 
age  by  the  probes  is  negligible.  Also,  at  a  Reynolds  number  of  11,000  a  dis¬ 
tance  of  1  cm  corresponds  to  y+  =  17,  permitting  high  spatial  resolution. 

The  use  of  oil  as  working  fluid  leads  to  higher  temporal  resolution  of  fluctu¬ 
ating  components  than  in  the  common  fluids,  air  and  water.  Velocity  measure¬ 
ments  are  with  commercial  hot  film  probes  which  are  calibrated  via  a  time-of- 
f light  technique  moving  the  support  cart  along  the  channel. 

For  measurements  with  laterally-converging  flow,  the  channel  floor  was 

o 

modified  to  a  convergence  angle  of  2  (Figure  5.1b)  giving  a  potential  inter- 

—  6  —6 

mediate  range  of  acceleration  parameters  1  x  10  <  <  3  x  10  as  desired 

to  examine  conditions  where  the  turbulence  structure  was  expected  to  be  modi¬ 
fied.  The  channel  was  straightened  and  reinforced  and  a  new  traversing  mech¬ 
anism  was  designed  and  constructed  for  measurements  in  the  entrance  half  of 
the  channel.  Data  acquisition  programs  were  modified  for  digital  application 
of  the  hot  film  calibration  as  a  function  of  the  oil  temperature  at  low  overheat 
ratios.  Experimental  studies  of  flow- induced  vibrations  were  conducted  in 
order  to  design  probe  supports  for  measurements  deep  in  the  channel. 

The  lateral  convergence  is  induced  by  a  ramp-shaped  insert  of  2°  which  has 
been  placed  in  the  bottom  of  the  open  channel,  as  shown  in  Figure  5.1b.  The  oil 
passes  along  a  return  channel  of  the  same  length  and  approximate  wall  spacing(s) 
giving  L/s  ~  30  for  flow  development  before  passing  through  a  series  of  curved 
vanes  in  the  180°  bend,  which  forms  the  entrance  to  the  test  channel.  With  a 
typical  depth  of  80  cm  the  aspect  ratio  is  about  3.6  for  the  open  channel;  with 
the  free  surface  approximating  a  zero  shear  surface,  this  configuration  is  com¬ 
parable  to  a  closed  channel  with  an  aspect  ratio  of  seven.  This  geometry  gives 
mildly  three-dimensional  velocity  distributions  away  from  the  edge  region 
(i.e. ,  bottom) . 

In  the  test  channel  the  entrance  region  is  Ax/s  ~  15  for  flow  development 
ahead  of  the  convergence  section  of  length  Ax/s  ~  17  to  the  minimum  cross 
section  ("throat").  The  main  measurements  in  the  converging  section  are  taken 
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at  x/s  =  31  (or  Ax  =  s  ahead  of  the  throat),  where  a  wall  sensor  is  also 
mounted . 

For  direct  comparison,  measurements  approximating  fully-developed  condi¬ 
tions  (Ky  -  0)  were  obtained  by  removing  the  ramp  from  the  bottom  of  the 
channel.  Otherwise,  locations  remained  the  same. 

6.  TECHNIQUES  AND  CALIBRATIONS 

The  techniques  and  calibrations  employed  were  essentially  those  of 
EcVelmann  [1974]  except  as  noted. 

6.1  Data  acquisition 

At  the  maximum  oil  velocity  of  23  cm/sec,  giving  Re  ^  11,000,  the  hot  film 
signals  are  sampled  for  over  one  hour  at  50  Hz  to  obtain  mean  values  within  one 
per  cent.  Maximum  significant  turbulence  frequencies  are  less  than  5  Hz  and 
bursting  frequencies  are  estimated  to  be  of  the  order  of  0.1  Hz  or  less. 

Time  series  data  were  acquired  via  a  DEC  PDP-15  system  from  four  hot-film 
sensors  simultaneously;  an  X-probe  for  u(t,y)  and  v(t,y),  a  wall  element  for 

r  (t)  and  a  single,  fixed  upstream  sensor  to  serve  as  a  reference.  Measure- 

w  + 

ments  were  taken  in  the  viscous  layer  (y  ~  5,  7,  10,  15,  25)  and  at  the  center- 

plane  for  Ky  =  1.6  x  10  6  and  Ky  =  2.4  x  10  6  for  the  cases  with  stream- 
wise  acceleration.  These  data  are  stored  on  digital  magnetic  tape  in  the  compact 
format  produced  by  the  VAX  "backup"  utility.  Some  additional  measurements  of 
mean  turbulence  quantities  for  fully-developed  flow  in  the  same  channel  are 
available  at  the  same  Reynolds  numbers  for  comparison  [Eckelmann,  1974; 

Randolph,  1983]. 

The  mean  signals  are  approximately  suppressed  by  a  constant  voltage  and 
then  are  amplified  to  improve  the  accuracy  (so-called  "buck  and  gain"  technique) 
before  analog-to- digital  conversion  and  storage.  Oil  temperature  is  measured 
manually  with  a  glass  thermometer  and  by  "cold”  resistance  measurements.  For 
the  X  probe  the  overheat  ratio  is  set  at  about  five  percent  (rather  than  the 
more  common  one  or  two  percent)  in  order  to  reduce  sensitivity  to  gradual 
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temperature  variations.  Despite  control  of  the  room  temperature,  these  pre¬ 
cautions  are  necessary  to  reduce  the  uncertainties  in  the  velocities  to  below 
one  percent. 

In  data  reduction,  a  subroutine  is  employed  to  apply  a  correction  for  the 
(small)  difference  between  operating  temperature  and  calibration  temperature, 
to  calculate  the  velocity  across  the  sensor  and,  for  the  X-probe,  to  apply  a 
calibrated  relation  for  angular  dependence  to  deduce  the  instantaneous  velocity 
components.  These  time  series  of  velocity  components  are  integrated  to  yield 
the  various  mean  statistics  desired. 

6.2  Velocity  calibration 

For  both  velocity  probes  the  sensors  were  calibrated  by  the  time-of-f light 
technique  of  Ecklemann  [1970,  1974).  The  support  cart  for  the  main  station  was 
linked  to  the  one  for  upstream  measurements  with  a  sturdy  bar  linkage,  so  the 
calibration  runs  were  conducted  simultaneously. 

Microswitches  1  meter  apart  were  connected  to  an  electronic  timer  to  measure 
the  time-of-f light  over  a  calibrated  distance.  Acceleration  to  a  steady  ve¬ 
locity  was  accomplished  before  the  cart  reached  the  first  switch  and  decelera¬ 
tion  occurred  after  the  last.  During  the  motion  of  the  cart  between  the  micro¬ 
switches  the  sensor  voltages  were  measured  with  integrating  digital  voltmeters. 

Oil  temperature  was  measured  with  a  glass  thermometer.  Although  it  took  two 

to  three  hours  to  obtain  the  desired  ten  to  fifteen  calibration  points,  the 

temperature  drift  during  this  period  was  negligible.  Therefore,  the  calibration 

temperature  was  taken  as  the  reference  temperature  for  the  resulting  velocity 

calibration  curves  V  , (e)  for  each  sensor. 

cal 

A  least  squares  approximation  procedure,  developed  by  Randolph  [1983],  was 
employed  to  determine  the  coefficients  ai  in  fourth-order  calibration  curves 
of  the  form 

2  3  4 

V  =  a  +  a  e  +  a  e  +  ae  +  a  e 
1  2  3  4  s 
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Weighting  of  the  calibration  points  was  adjusted  to  yield  agreement  within  1/2 
percent  or  better  for  the  difference  between  the  resulting  curve  and  any 
measured  point. 

6.3  Temperature  correction 

During  the  measurements  with  flow,  the  oil  temperature  was  close  to  the  cali¬ 
bration  temperature  but  not  necessarily  the  same,  so  a  correction  was  applied  for 
the  slight  difference.  Development  of  the  correction  factor  is  summarized  in 
section  7.3  concerning  the  wall  shear  sensors. 

The  oil  temperature  and/or  equivalent  cold  resistance  of  the  sensor  was 
measured  before  and  after  the  data  acquisition  for  the  run.  For  the  correction 
the  instantaneous  measured  voltage  was  adjusted  to  the  value  it  would  have  been 
at  the  same  velocity  if  the  oil  had  been  at  the  reference  temperature.  This 
first  order  correction  was  of  the  form 

e  R  R  r 

ref  1  c.op  -  c.ref 

e  1  +  2  R  -  R 

op  op  c ,  op 

Where  R£  stands  for  cold  resistance  and  the  subscript  op  denotes  the  measuring 
or  operating  conditions  of  the  actual  run.  The  correction  is  actually  intro¬ 
duced  by  adjusting  the  values  of  the  gain  and  bucking  voltage  which  the  data 
reduction  programs  employ  to  convert  the  integer  digital  sample  records  into 
equivalent  (floating  point)  voltages. 

6.4  Velocity  calculation 

In  some  preliminary  measurements  a  single  sensor  probe  was  used  to  measure 
the  streamwisa  component  u.  In  that  case,  the  velocity  was  calculated  directly 
from  its  calibration  relation  V  ^Ce)  after  measuring  the  anemometer  volt¬ 
age.  For  the  X-probe  a  cosine  approximation  was  applied  via  an  algorithm 
employed  by  Randolph  [1983], 

u  =  V  ,  cos  a  +  V  cos  a 
cal,i  i  cal, 2  2 

v  =  (V  -  V  ,  )  cos  a  cos  a 

cal , i  cal ,  2  i  2 
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where  Vcg^  ^  represents  "calibration  velocity"  after  correcting  for  temp¬ 
erature,  converting  to  voltage  and  applying  the  calibration  function  V^^Ce). 

The  angles  of  the  two  wires  were  measured  to  be  +  45.0  degrees  by  Herr  Nortemann 
who  constructed  and  repaired  hot  film  sensors  as  appropriate. 

In  the  present  work  the  normal  mean  velocity  v  was  not  adjusted  to  zero, 
whereas  earlier  workers  with  fully  developed  flows  did  "correct"  the  sensor 
angles  to  yield  v  =  0  before  further  processing  [Randolph,  1983;  Brodkey,  1990). 
We  retained  the  non-zero  result  since  v  is  not  necessarily  zero  in  flows  with 
streamwise  acceleration.  The  deduced  fluctuation  v  is  taken  around  v  in  the 
same  manner  that  the  streamwise  fluctuation  is  determined.  Since  v  is  cal¬ 
culated  as  a  small  difference  between  larger  quantities,  its  percent  un¬ 
certainty  is  considerably  more  than  that  of  u.  The  non-zero  values  of  v  cal¬ 
culated  for  the  fully  developed  flows  (i.e.,  without  ramp)  are  consequences  of 
that  uncertainty  and  an  apparent  wall  interaction  observed  by  Randolph  [1983]. 
Therefore,  for  the  accelerating  flows,  in  order  to  decide  whether  non-zero 
values  of  v  are  meaningful,  it  is  necessary  to  compare  to  the  fully-developed 
results . 

6.5  Wall  distance,  y 

The  theodolite  technique  of  Eckelmann  [1970,  1974]  was  employed  to  set  the 
reference  location  of  the  crossing  of  the  sensors  of  the  X-probe  to  within 
about  +  0.1mm.  For  y  <  30mm,  a  dial  indicator  with  an  estimated  reading 
uncertainty  of  +  0.002mm  was  used  to  determine  transverse  locations  relative 
to  this  reference.  Beyond  that  distance  a  vernier  rule  was  employed;  it  could 
be  read  to  ±  0.1mm. 

The  wall  distance  y  (or  y+)  is  typically  used  as  a  coordinate  on  the  profile 
plots  and  to  identify  approximately  equivalent  locations  when  comparing  acceler¬ 
ated  flows  with  fully  developed  flows.  For  these  purposes  high  accuracy  is  not 
necessary.  While  Eckelmann  [1974]  determined  the  location  y  =  0  by  linear  ex¬ 
trapolation  of  his  near  wall  data,  the  present  report  uses  physical  measurement 
of  the  length  of  the  pin  which  served  for  the  reference  location  in  the  theo¬ 
dolite  measurements.  A  slight  correction  for  displacement  of  the  wall  surface 
(due  to  pressure  differential  between  the  two  parallel  flow  channels.  Fig.  5.1) 
is  also  included. 


-20- 


Precision  is  required  for  the  location  of  the  one-point  velocity  measurement 
IMeKligot,  1984]  applied  to  determine  mean  wall  shear  stress  (to  calculate  y+ 
and  to  calibrate  the  wall  sensor).  Sensitivity  studies  showed  a  typical  influ¬ 
ence  coefficient  of  the  order  of  unity  (see  Figure  7.2  later),  i.e.  a  one  per¬ 
cent  error  in  y+  would  cause  a  one  percent  error  in  the  value  of  t  deduced. 

+  w 

For  the  worst  case  in  this  determination,  the  uncertainty  in  y  was  about  0.04 
percent  so  the  uncertainty  in  the  velocity  measurement  dominates  the  wall  shear 
stress  determination. 

6.6  Measurement  of  pressure  gradient 

With  oil,  the  pressure  gradient  corresponding  to  a  significant  value  of  K  is 

P 

miniscule.  Thus,  it  was  difficult  to  obtain  an  accurate  measurement  of  dp/dx  for 
determining  and  K^.  Several  methods  were  considered  and/or  attempted. 

Static  pressure  taps  were  designed  and  installed  by  Reichardt  along  the  wall 
in  the  original  channel.  These  were  no  longer  useful. 

Two  pitot-static  probes  were  installed  in  the  channel  approximately  in  the 
same  manner  as  the  X-probe,  but  displaced  longitudinally.  The  difference  be¬ 
tween  the  static  pressure  taps  on  the  two  probes  was  measured  with  an  MKS 
Baratron  Differential  Pressure  Meter.  We  used  the  most  sensitive  transducer, 
one  with  a  span  of  1  mm  Hg.  Unfortunately,  the  differences  in  the  sizes  of  the 
transducer  plenums,  the  length  of  the  necessary  tubing  (volume  considerably 
greater  than  transducer)  and  slight  variation  in  the  room  temperature  led  to 
slow  oscillations  of  the  measured  differential  pressure.  The  amplitude  of 
these  oscillations  was  greater  than  the  magnitude  of  the  pressure  drop  so  the 
approach  proved  to  be  impractical. 

Ultimately  a  hooR  gage  technique  was  applied.  Essentially,  the  drop  in  the 
oil  surface  elevation  due  to  flow  was  measured  at  several  locations  for  the 
flow  rates  of  interest.  The  pressure  gradient  dp/dx  was  then  approximated  as 
Az/Ax  with  flow.  Precision  of  the  individual  measurements  was  about  0.05 
mm  which  propagated  to  an  estimated  experimental  uncertainty  of  0.14  mm  for  Az 
(the  per  cent  uncertainty  of  Ax  was  negligible).  Since  measurements  of  Az 
varied  from  0.15  to  0.5  mm  for  the  range  of  interest,  the  estimated  experi- 
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mental  uncertainty  in  dp/dx  and  is  large  (100  to  30  per  cent, 
respectively! ) . 

The  consequent  estimates  of  K  can  be  considered  to  be  little  better  than 

P 

order-of-magnitude  estimates.  However,  their  consistency  gives  confidence 
that  meaningful  comparisons  of  the  various  results  can  be  made  in  order  to 
examine  the  effects  of  significant  increases  in  the  (non-dimensional)  stream- 

wise  pressure  gradient.  Fortunately,  the  determination  of  t  is  insensitive 

w 

to  the  estimated  experimental  uncertainty  in  dp/dx. 

6.7  Amplifier  Rain 

The  amplifier  gain  was  calibrated  in  situ.  That  is,  it  was  determined  for 
the  same  circuitry  and  settings  as  were  used  for  the  actual  measurement  runs. 
With  the  instrumentation  ready  for  a  run,  the  leads  from  the  anemometers  were 
disconnected  and  replaced  by  calibrated  voltage  sources.  The  digital  data 
acquisition  program  was  then  run  as  if  it  were  a  normal  run  and  the  amplifi- 
tion  factors  were  calculated  from  the  mean  signals  recorded. 

The  same  procedure  was  employed  to  deduce  the  slight  zero  error  that  remained 
after  the  initial  adjustments  of  the  amplifiers.  The  inputs  to  the  amplifiers 
were  shorted  and  the  normal  digital  data  acquisition  plus  integration  of  the 
recorded  signals  then  yielded  the  amplifier  zero  error.  Typically,  it  was  less 
than  one  or  two  digits  in  1024,  but  it  was  accounted  for  in  the  data  reduction. 

A  short  program  (DINP.FOR)  calculated  the  temperature  correction  and  then 
the  input  values  of  gain  and  bucking  voltage  used  by  the  main  data  reduction 
programs,  DUVTR . FOR ,  DCORR.FOR  and  DSPEC.FOR. 


6 . 8  Channel  width 


The  channel  width  was  measured  along  its  length  by  Herr  Nortemann  with  a 
large  dial  indicator.  The  estimated  experimental  uncertainty  was  +  0.05  mm. 
Before  the  experiment  was  conducted  the  spacers  above  the  oil  surface  were  reset 
to  give  equal  spacing  along  the  channel.  At  the  measuring  position,  the  deflec- 


-22- 


tion  near  the  oil  surface  was  approximately  1.5  mm  from  no  flow  to  maximum  flow. 
This  effect  was  included  in  the  calculation  of  the  probe  location  y  and  its 
uncertainty. 


The  orientation  of  the  channel  wall  relative  to  the  track  of  the  calibration 
cart  was  also  checked  by  mounting  a  dial  indicator  to  the  probe  support  struc¬ 
ture  (Figure  5.1)  and  moving  it  along  the  channel. 

7.  WALL  SHEAR  STRESS 

The  one-point  technique  suggested  by  HcEligot  [1984;  1985]  was  extended 
slightly  and  was  applied  to  determine  the  wall  shear  stress.  The  wall  sensor 
calibration  was  done  essentially  in  situ.  With  the  X-probe  at  y+  =  10  it  was 
away  from  significant  wall  interaction  but  close  enough  that  the  uncertainty  due 
to  choice  of  a  turbulence  model  was  minimal. 

Typically  a  non-linear  polymonial  was  employed  to  relate  the  instantaneous 
wall  shear  stress  to  the  instantaneous  voltage  signal 

2 

t  =  a  +  a  e  +  a  e  +  . 

w  1  2  3 

at  the  temperature  ("cold  resistance")  chosen  for  the  reference  calibration. 

The  slight  difference  between  the  reference  temperature  and  temperature  measured 
during  a  run  was  corrected  when  the  calibration  was  applied  to  deduce  T^Ct). 


7.1  Determination  of  t  for  calibration 
_ w _ 

The  one-point  technique  deduces  the  mean  wall  shear  stress  from  the 
measured  mean  streamwise  velocity  u  at  a  measured  distance  from  the  wall, 
via  a  turbulence  model  for  the  viscous  layer  (Figure  7.1).  McEligot  [1984] 
presented  the  details  including  sensitivity  to  turbulence  model  and  measure¬ 
ments  (Figure  7.2),  resulting  experimental  uncertainties  (Figure  7.3),  and 
the  Couette  flow  analysis  upon  which  it  was  based.  For  fully  developed  flows 
the  assumption  that  the  convective  terms  are  negligible  in  the  governing 
momentum  equation, 
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u  3u  v  3u 
3x  +  3y 


1  d£  1  3t 
p  dx  p  3y 


is  ideally  approached.  However,  for  spatially  accelerating  flows  the  relative 
magnitude  of  these  terms  depends  on  the  applied  pressure  gradient  and  the 
distance  y  from  the  wall  [Hanratty,  1988]. 

To  account  for  the  convective  terms,  we  adopted  a  suggestion  of  Finnicum  and 
Hanratty  [1988]  based  on  earlier  work  of  Julien  et  al.  [1969].  For  Couette 
flow  the  local  shear  stress  -r(y)  is  given  as 


Ay)  -=  1  +•  K  y+ 
P 


'f'  3 

where  t  (y)  =  t.  . /t  here  and  K  =  (v/pu  )dp/dx.  Given  a  reasonable 
(y)  w  p  t 

turbulence  model,  the  one-point  technique  effectively  determines  the  value  of 

t  which  satisfies 

w 


u  ,f 


”  ~  Jo  -efffy) 


where  the  subscript  m  refers  to  the  measuring  point  and  represents  the 


turbulence  model. 


For  sink  flows,  Julien  et  al.  [1969]  developed  the  relation 


c  y 

+ .  +  I  -  f  1 _  /*  ,  +x2  *  + 

t  <y  )  -  1  +  K  y  1  -  —  7  (u  )  dy 


to  account  for  the  convective  terms.  We  implemented  this  revision  by  multi¬ 
plying  the  term  K  y+  in  existing  program  TAU.FOR  [McEligot,  1984]  by  the 
grouping  in  brackets  above.  For  fully-developed  flow  this  change  is  sup¬ 
pressed.  The  integral  which  makes  the  procedure  implicit  was  evaluated  by 
successive  substitution  during  the  iterative  process  for  or  u^.  The 
revised  program,  TAU2.F0R,  is  listed  in  Appendix  C. 
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The  effect  of  the  convective  terms  was  evaluated  for  several  accelerated 
flows  of  interest  by  conducting  the  calculations  with  and  without  the  correc¬ 
tion.  Typical  results  were 


Kp 

=  0.011 

=  0.011 

=  0.02 

=  0.02 

y+ 

10 

25 

10 

15 

Aut 

0.1% 

0.8% 

0.5% 

0.9% 

AK 

0.3% 

2.5% 

1.5% 

2.6% 

P 


Since  values  with  the  X-probe  at  y+  =  10  are  used  to  determine  for  calibra¬ 
tion  and  for  calculating  wall  coordinates,  the  inclusion  of  the  convective 
terms  led  to  no  great  change. 


7.2  In  situ  calibration 


As  explained  in  the  section  above,  during  calibration  of  the  wall  sensor 
the  mean  wall  shear  stress  is  determined  by  the  one-point  method  [McEligot, 

1984]  as  revised.  Actually,  the  calibration  data  are  extracted  from  normal 
measuring  runs.  Thus,  the  flow  is  turbulent  with  the  same  fluctuating  field 
present  during  calibration  as  during  application. 

The  calibration  relation  for  the  fluctuating  wall  shear  stress  is  taken  to  be 

t  (t)  =  a  +•  a  e  +  a  e2  +  .  .  .  . 

W  12  3 

where  e  is  the  instantaneous  voltage  signal  from  the  sensor.  By  integration 
with  respect  to  time,  one  finds  the  relation  for  mean  values  to  be 

-  ~  ~7 

t  =  a  +  a  e  +-  a  e  +  .... 
w  1  2  3 

The  data  reduction  program  for  mean  velocities,  DUVTR,  was  modified  slightly  so 
that  with  suitable  choices  of  calibration  coefficients  it  would  calculate  the 

mean  voltage  statistics  of  the  wall  sensor,  e,  e7,  e7  and  e7  as  part  of  the 
data  reduction  for  a  normal  run. 

In  order  to  calculate  the  calibration  coefficients  a^  for  an  n-th  order 
polynomial,  one  needs  calibration  data  at  n  t  1  (or  more)  levels  of  t  .  The 
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values  are  then  determined  by  solving  the  n  +  1  simultaneous  equations  for 

t  .  or  by  a  least-squares  fit. 
w  #  i 

7.3  Temperature  correction 


Although  the  temperature  of  the  room  was  closely  regulated  by  thermostatic 
control,  there  were  slight  variations  from  run  to  run  during  the  overall 
period  of  the  measuring  program,  say  +  0.1  C.  During  an  individual  run  the 
fluctuation  was  a  few  hundredths  of  a  degree  at  most.  Since  the  overheat 
ratio  used  for  the  wall  sensor  operation  was  relatively  low  (equivalent  to 
about  10  C) ,  a  first  order  temperature  correction  was  developed.  The  approach 
is  comparable  to  one  suggested  by  S.J.  Kline  a  number  of  years  ago  [reference 
citation  now  forgotten]. 

We  choose  a  standard  or  calibration  temperature  to  which  the  data  will  be 
referred.  The  calibration  relation  Tw(e)  is  then  developed  for  that  temper¬ 
ature.  In  converting  calibration  data  to  this  temperature  and  in  applying  the 
calibration,  the  key  question  is  —  for  the  shear  stress  existing  at  the  oper¬ 
ating  temperature,  what  voltage  signal  would  the  same  shear  stress  cause  at  the 
standard  reference  temperature? 


Figure  7.4  describes  the  conceptual  process  graphically.  The  individual 

curves  represent  the  functional  behavior  u(e)  expected  for  different  fluid 

temperatures  (or  temperature  differences  or  overheat  ratios)  while  the  sensor 

operating  resistance  is  held  constant.  The  variable  u  is  the  wall  shear  stress 

or  the  velocity,  as  appropriate,  and  e  is  the  required  sensor  voltage.  For 

example,  T  and  T  may  be  considered  the  fluid  temperatures  at  which  measure- 
1  2 

ments  are  made  while  T  is  the  temperature  chosen  for  the  reference  condition. 

3 

If  we  have  data  at  T  for  a  given  value  of  u  (and  e  ) ,  we  wish  to  determine 

i  i 

the  value  e  which  would  be  required  for  the  same  value  of  u  at  T  .  A  Taylor 

3  3 

series  approximation  would  give 


e  =  e 

3  1 


AT  + 
u 
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or 


&e  +. . . . 

u 

where  0  represents  a  temperature  difference,  overheat,  or  other  convenient 
function  of  temperature.  In  practise,  we  ultimately  form  the  ratio  e^/e^  80 
the  correction  is  a  multiplier  for  the  measured  voltage. 

The  first  order  correction  takes  the  steady  state  energy  balance  of  the 
sensor  to  be 


2 

r —  =  q  =  hA  (T 
R  conv  s  op 

op 


T  ) 

CD 


and  attributes  the  primary  temperature  effect  to  the  term  (T  -  T  )  and 
neglects  the  variation  of  other  quantities  due  to  temperature. 

The  resistance  R  is  held  constant  in  "constant  temperature"  operation 
°P 

of  the  anemometer.  With  the  temperature  coefiicient  of  resistance,  one 
can  write 


R  =  R  [1  +  a  (T  -  T  ) ]  =  R  [1+a  (T  -  T  ) ] 
op  o  op  o  c  op  00 

where  the  subscript  o  represents  a  reference  condition.  In  our  case  we  let 


R  =  R  =  R(T  )  =  R(T  ) 
o  c  o  00 


where  R  is  called  the  cold  resistance, 
c 


The  energy  balance  can  then  be  rearranged  to 


hA  R  ,  _  _  . 

2  s  op  (R  -  R  ) 
*  -  °r  C 
c 


-27- 


or,  with  a  definition  for  0, 


R  -  R 

e2  =  fn{u,  constants}  °PR  h  fn{u}  0 

c 


Further  rearrangement  gives 


e  =  /FnTuT  I©] 


ft  -  \  ^  ™'/2 


For  the  first  order  approximation 


,  1  A0  ,  ,  1  A0 

1  +  30  e  _1  20 

l  i 


e  ,  R  R  -  R 

_1  „  .  I  _2£  .  -£1 - £3 

e  2  R  R  -  R  , 

1  C  3  OP  Cl 

For  the  correction  of  velocity  signals  to  the  chosen  calibration  reference 
temperature,  we  approximated  R  /RCa  -  1  (along  with  Kline);  this  introduced  a 
small  percent  error  in  a  small  percent  correction  (typically  less  than  one 
percent).  For  the  wall  sensor  calibration  Rop^C3  was  included. 

For  more  convenient  application,  the  correction  may  be  recast  into  terms 
of  temperatures  via  the  relation 


R  as  R  +  —  (T  -  T  ) 
Cl  C  3  3T  1  3 
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The  result  is 


e 

3 

e 

i 


1 


1  -  (R  /R  )(1  +  B) 

C  3  Op 


where 


B 


CT 

i 


T  ) 

3 


1 

R 

C  3 


3R 
_ c 

3T 


7.4  Application  of  wall  shear  stress  calibration 

For  the  fully-developed  conditions,  measurements  are  available  at  four  wall 
shear  stress  levels  corresponding  to  runs  at  four  flow  rates  Cor  pump  speeds). 
A  temperature  of  24.8  C  was  chosen  to  serve  for  the  calibration  curve;  this 
temperature  was  in  the  midst  of  the  range  of  operating  temperatures.  The  mean 

—  — 7  — 3T  — y 

voltages  e,  e  ,  e  and  e  for  each  run  were  converted  to  the  values  that  would 
have  been  at  24.8  C  via  the  approach  of  the  previous  section. 

Due  to  the  limited  range  of  values  available  for  calibration  purposes, 
a  second  order  polynomial  was  used  for  the  actual  curve, 

2 

t  =  a  +  a  e  +  a  e 
W  1  2  3 

The  coefficients  were  determined  by  fitting  the  relations 

—  —  2 
T.=a  +  a  e.  +  a  e. 

W1  1  2  1  3  1 


through  the  data  for  u^  =  0.52,  0.80  and  1.04  cm/sec,  spanning  the  range  of 
measurements  to  be  reported. 

The  mean  t  data  and  the  resulting  calibration  curve  t  (e)  are  plotted  on 
w  w 

Figure  7.5.  It  is  seen  that  i^Ce)  is  near  linear  over  this  range.  Further, 

it  does  not  differ  significantly  from  Tw(e),  although  the  second  order  term 

2  ,  , 
ase  is  not  negligible. 

For  the  calculations  of  mean  statistics,  correlations  and  spectra  involving 

t  for  runs  at  the  various  X-probe  locations,  y+,  the  temperature  correction  to 
w 

24.8  C  was  introduced  by  correcting  the  input  gain  and  bucking  voltage  to  account 
for  e  /e  .  Then  the  coefficients  a  ,  a  and  a  (with  a  =  a  =0)  were  employed 

3  1  1  2  3  4  5 
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in  the  data  reduction  program  to  convert  each  voltage  sample  for  the  wall  sensor 
to  its  equivalent  instantaneous  wall  shear  stress. 

For  normalization  of  wall  shear  stress  quantities  in  correlations,  spectra, 

skewness,  flatness,  etc.,  the  mean  value  was  determined  by  integration 

of  the  instantaneous  values  of  tw  in  the  data  reduction  program.  However, 

wall  coordinates  such  as  y+  and  t+  were  calculated  from  t  (or  u  ) 

w  T 

deduced  from  the  one-point  velocity  technique  of  McEligot  [1984].  The  differ¬ 
ences  in  involved  were  slight  (and  at  y+  -  10  essentially  nil  since 
it  was  the  basis  of  the  calibration). 

7.5  Treatment  of  wall  shear  stress  signal  for  accelerated  flows 

The  same  wall  sensor  was  used  in  the  same  position  in  the  laterally  con¬ 
verging  configuration  as  in  the  fully  developed  version.  Therefore,  the  basic 
thermal  behavior  of  the  sensor  should  be  the  same  for  the  two  cases  and  the 
calibration  determined  for  fully  developed  flow  should  be  the  same  as  for  the 
accelerated  flows.  However,  the  operating  resistance  of  the  sensor  and  its 
cold  resistance  were  slightly  different  in  the  two  cases.  Since,  for  the 
accelerated  case,  complete  details  of  the  gain,  bucking  voltage  and  cold  resist¬ 
ance  at  the  measuring  temperature  are  now  available  for  only  one  strain  rate  or 
wall  shear  stress  level  (thanks  to  Lufthansa  baggage  handling  in  Frankfurt),  it 
is  necessary  to  combine  these  data  with  the  calibration  from  the  fully 
developed  runs. 

For  steady  or  quasi-steady  state,  the  energy  balance  for  the  wall  sensor 
can  be  written  as 

e2/R  =  hA  (T  -  T  ) 
s  w  c 

where  the  heat  transfer  coefficient  represents  heat  transfer  to  the  fluid  by 
three  modes, 

h  =  +  h  .  +  hc 

NC  eff,  cond  fc,t 

buoyancy,  conduction  via  the  substrate  and  forced  flow  effects  which  vary  with 
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the  wall  shear  stress.  It  is  the  variation  of  the  last  which  provides  the 
usual  calibration  of  voltage  signal  versus  wall  shear  stress.  If  the  Leveque 
solution  [1928]  describes  the  forced  flow  contribution  and  the  temperatures  (or 
resistances)  are  held  constant,  the  relation  can  be  rearranged  to 


1/3 


w 


c  e  -  c  h„_  -chrc  .=ce  -A 
1  2  NC  3  eff.cond  1 


with  the  c^  representing  constant  coefficients.  Over  the  range  of  the 
present  experiments  this  relation  may  be  approximated  as 


2 

t  =  a  +  a  e  +  a  e 
w  1  2  3 


as  shown  in  section  7.4  above. 


It  is  expected  that,  aside  from  the  change  in  temperature  difference  or 

overheat  ratio  (which  is  treated  as  in  section  7.3)  the  main  effect  of  a  change 

in  sensor  operating  resistance  R  is  to  modify  the  "constant"  contributions  from 

op 

buoyancy  effects  and  substrate  conduction.  For  application  to  the  accelerated 
data,  the  values  of  and  a ^  are  taken  from  the  three  point  calibration 
(section  7.4)  and  coefficient  a^  is  calibrated  via  the  data  with  the  probe 
at  y+  -  10  and  the  lower  flow  rate  (corresponding  to  a  pump  speed  of  1200  rpm) . 

The  time  series  for  the  data  at  this  point  are  integrated  to  give  the  mean, 
voltage  quantities  for  the  relation 


t  =  a  +  a  e  -  a  e 
w  1  2  3 


and  the  mean  wall  shear  stress  is  determined  from  the  mean  velocity  u(y+  =  10) 
in  accordance  with  the  procedure  of  section  7.1.  Thus,  the  desired  coef¬ 
ficient  is  given  as 


a  =  t  -  ae-ae 
1  w  2  3 


from  the  calibration  measurements. 
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For  the  laterally  converging  flows,  the  temporal  wall  shear  stress  is  normal¬ 
ized  by  its  mean  value  for  application,  thereby  reducing  the  effects  of  un¬ 
certainties  in  this  calibration  procedure.  For  wall  coordinate  scaling  as  in 

determining  u+,  y+,  etc.,  the  one  point  velocity  method  (section  7.1)  is  used  to 

obtain  t  . 

w 

When  the  coefficients  a  and  a  from  from  the  fully-developed  calibration 

2  3 

were  used  to  deduce  a  ,  it  was  found  that  the  resulting  calculated  value  of  t 

i  w 

was  reasonable.  However,  the  rms  fluctuation  of  t  appeared  to  be  much  too 

w 

small  (about  1.2  percent  of  t  ) .  This  difficulty  was  traced  to  the  second 

w 

order  polynomial  in  the  calibration.  While  the  variation  of  was  within  its 
calibrated  range,  the  voltages  were  lower  than  those  in  the  fully- developed  cali 
bration.  The  values  of  a  and  a  led  to  a  minimum  in  the  calibration  curve 

2  3 

t  (e)  near  the  mean  value  e  for  the  accelerated  flow.  Thus,  significant  fluctua 

w 

tions  in  e  yielded  only  slight  variations  in  the  calculated  t  calculated 
t  w 

and  a  low  rms  value  of  t  was  deduced. 

w 

To  counter  this  difficulty  a  linear  calibration 

TW  =  ^2® 

was  ultimately  employed.  The  coefficient  a  was  deduced  from  the  slope  of  the 

2 

fully-developed  correlation  between  its  lowest  and  highest  values  of  t  .  As 

noted  earlier,  the  relation  t  (e)  was  approximately  linear  over  this  range 

w  .  .  + 

(Figure  7.5).  The  constant  a^  was  then  adjusted  to  fit  the  data  at  y  =10  for 

the  measurements  at  K  =  -0.02  as 

P 


a  =  t  -  a  e 
1  w  2 


The  deduced  value  of  rms  then  was  approximately  twenty  percent  of  and  the 

skewness  and  flatness  factors  S  and  F  showed  magnitudes  comparable  to  S 

t  t  u 

and  F  . 
u 


7.6  Determination  of  pressure  gradient  parameter. 


K 

_E 


The  pressure  gradient  parameter  is  defined  as 
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v  dp 
p  pu  dx 


The  friction  velocity  u  was  determined  via  program  TAU2.F0R  along  with 
as  described  above.  The  pressure  gradient  dp/dx  was  estimated  with  the  "hook 
gauge"  measurements  as  mentioned  in  section  6.6  and  its  experimental  uncertainty 
is  relatively  large.  However,  even  though  dp/dx  is  an  input  in  the  determina¬ 
tion  of  u  ,  the  resulting  value  of  u  is  insensitive  to  the  uncertainty  of  dp/dx 

T  T 

or  K  (Figure  7.2).  Consequently,  the  estimated  experimental  uncertainty  of 
P 

K  is  approximately  that  of  dp/dx. 

p 


8.  RESULTS 

The  primary  experimental  results  are  the  simultaneous  measurements  with  the 

wall  sensor  for  t  and  X-probe  for  u  and  v.  Velocity  data  were  concentrated 
w 

in  the  viscous  sublayer  at  y  -  5,  7,  10,  15  and  25  plus  the  centerplane. 
Typically,  200,000  samples  of  each  signal  were  acquired  for  the  time  series. 
Four  sets  of  operating  conditions  were  selected  as  follows: 


Ramp  .  us 

.  .  K  =  — j —  Re  =  — 

Angle  p  pu^  dx  s  v 

"Fully-developed"  flows 
0  -  0.0083  7400 

0  -  0.011  5600 


„  v  dV 

v  vr  dx 

D 


0 

0 


Streamwise  accelerating  flows 
2°  -  0,011  8300  1.6  x  10 

2°  -  0.02  5600  2.4  x  10 


This  selection  of  conditions  permits  separating  effects  of  K  from  Re  and  from 

P  s 

K^.  For  example,  at  one  (approximately  constant)  Reynolds  number  there  are 
data  at  two  different  values  of  K  ,  and  vice  versa. 


Mean  turbulence  statistics  are  tabulated  in  Appendix  A.  Additional  mean 

velocity  profile  measurements  are  listed  in  Appendix  B;  these  are  required  for 

the  estimation  of  V,  in  the  calculation  of  K  .  Results  calculated  to  date  are 

b  v 
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these  mean  statistics,  selected  auto-  and  cross- correlations  plus  preliminary 
power  spectra  of  the  signals.  This  chapter  presents  those  results  with  dis¬ 
cussion  as  appropriate.  Usually  our  presentation  is  in  terms  of  wall  co¬ 
ordinates  . 

8.1  Mean  velocity  profiles 

The  streamwise  mean  velocity  profiles  are  presented  in  Figure  8.1  in  terms  of 

wall  coordinates.  These  data  confirm  the  trends  predicted  in  Figure  7.1;  as  K 

+  +  ^ 
increases  in  magnitude  the  curves  u  (y  )  also  have  highei.  values  away  from  the 

wall.  This  is  one  effect  of  streamwise  pressure  gradient  on  mean  turbulence 
structure.  It  corresponds  to  a  thickening  of  the  layer  dominated  by  viscous 
effects  as  in  the  mixing  length  wall  models  of  McEligot,  Ormand  and  Perkins 
[1966],  Huffman  and  Bradshaw  [1972],  Launder  and  Jones  [1969]  and  others.  The 
trend  is  also  predicted  by  the  direct  numerical  simulations  of  Spalart  [1986]. 
Comparable  data  have  been  obtained  by  Jones  and  Launder  [1972]  and  Loyd,  Moffat 
and  Kays  [1970]  for  transversely  converging  flows,  by  Perkins  and  McEligot 
[1975]  and  Shehata  [1984]  for  strongly  heated  gas  flows  in  tubes  and  in  turbu¬ 
lent  flows  with  drag  reduction  by  polymer  additives  [Berman,  1977;  Harder  and 
Tiederman,  1989]. 

The  deduced  values  of  the  transverse  mean  velocity  v+  are  tabulated  in 
Appendix  A  for  the  viscous  layer  and  the  centerplane.  These  data  can  be  sub¬ 
ject  to  a  variety  of  difficulties.  The  calculation  of  v  from  the  X-probe 
signals  involves  determining  the  small  difference  between  two  much  larger 
signals.  (For  a  fully-developed  flow,  v  is  ideally  zero  and  consequently 
percent  experimental  uncertainties  are  infinite.)  In  a  fully  developed  flow  in 
the  same  oil  channel  Randolph  [1988]  found  that  as  his  X-probe  approached  the 
wall,  the  measurements  yielded  flow  angles  which  were  not  parallel  to  the  wall 
—  a  possible  indication  of  wall-probe  interaction.  Moin  and  Spalart  [1988] 
modeled  the  behavior  of  a  typical  X-probe  via  direct  numerical  simulations  of  a 
turbulent  flow;  they  found  that  neglecting  the  spanwise  component  w  led  to 
larger  errors  in  v  than  in  u  and  the  separation  of  the  two  films  can  have 
strong  effects.  Moi.i  and  Spalart  did  not  present  calculations  of  the  expected 
measured  value  for  the  mean  v*(y+),  but  in  general  the  predicted  errors  in- 
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creased  as  y+  approached  zero  with  dependence  on  w  and  the  calibration  rela¬ 
tions  employed. 


Thus,  the  statistics  involving  v  are  best  considered  in  a  comparative  or 
relative  manner  rather  than  placing  too  much  confidence  in  their  absolute 
values.  That  is,  we  recommend  examining  the  effects  of  by  using  the  data 
for  the  lowest  value  as  a  reference  (fully-developed  flow,  =  -0.008) 

and  then  comparing  the  data  for  other  Kp  to  it.  Other  investigators  have 

arbitrarily  adjusted  their  sensor  angles  to  force  v  =  0  in  order  to  avoid  the 

appearance  of  this  type  of  problem.  As  mentioned  earlier,  for  a  laterally 
converging  flow  v  is  not  identically  zero  so  we  carry  the  values  actually 

calculated  (and  look  for  relative  effects  at  the  same  locations). 


The  mean  values  v(y  )  are  all  towards  the  wall  except  centerplane  values  for 
the  fully-developed  flows.  The  largest  values  appear  at  y+  *  7  to  10  and  then 
the  magnitudes  decrease  as  y+  increases.  However,  they  are  generally  less  than 
about  five  percent  of  the  value  of  u  at  the  same  y+.  Although  there  are  differ¬ 
ences  between  the  values  of  v*(y+)  for  the  various  runs,  there  is  no  clear  trend 
as  K  varies  and  the  difference  at  a  given  y+  is  generally  within  the  experi- 
ental  uncertainty.  The  major  variation  observed  is  with  y  and  is  likely  to  be 
a  consequence  of  the  difficulties  discussed  above. 


8.2  Root  mean  square  fluctuations 

The  normalized  values  of  the  r.m.s.  velocity  components  are  plotted  in 
Figures  8.2  and  8.3.  For  u'  the  maximum  value  approaches  three  at  y+  near  15 
in  agreement  with  the  data  of  Eckelmann  [1974]  and  others.  Also,  v'  is  much 
less  than  u’  and  continues  to  increase  within  the  viscous  layer  rather  peaking 
there  as  u'  does.  The  magnitude  of  v’(y+)  is  approximately  the  same  as  the  data 
of  Eckelmann  but  does  fall  slightly  lower. 

Eckelmann  compared  the  data  at  two  Reynolds  numbers,  corresponding  to  our  two 
fully-developed  conditions.  When  non-dimensionalized  with  wall  coordinates 
(his  Fig.  13)  the  data  collapsed,  i.e.,  there  was  no  apparent  effect  of  Reynolds 
number.  Our  Figures  8.2  and  8.3,  using  comparable  wall  scaling,  also  show  no 
significant  trends  with  variation  of  Kp.  That  is,  for  these  measures  of  the 
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turbulence  structure  there  appears  to  be  no  important  effect  of  the  stream- 
wise  pressure  gradient. 

8.3  Reynolds  shear  stress 

With  the  dimensions  of  the  channel  and  the  flow  rates  yielding  centerline 
values  of  y+  in  the  range  of  140  to  210,  the  total  shear  stress  varies  observ¬ 
ably  across  the  viscous  layer.  Therefore,  the  popular  constant  shear  layer 
approximation  is  not  valid,  even  for  the  fully-developed  flows. 

Near  the  wall  the  total  shear  stress  variation  is  given  approximately  as 


-r(y)  _  3u*  _  uv 

t  +  2 

w  dy  u 

T 


1  +  K  y+ 
P 


neglecting  the  convective  contribution  in  the  x-momentum  equation.  For  fully- 

developed  flow  between  parallel  plates  this  relation  is  exact  while  for  accel- 

+ 

erating  flows  it  will  diverge  as  y  increases  (see  the  earlier  discussion  in 
section  7.1  for  further  details). 


At  the  wall  u  =  v  =  0  so  uv  is  likewise  zero  there.  Since 


f(y) 


puv. 


the  Reynolds  shear  stress  increases  with  y  approaching  the  total  shear  stress 
relation  as  the  turbulent  transport  becomes  dominant  relative  to  viscous 
effects. 


Figure  8.4  presents  the  approximate  asymptotic  values  for  T(y  )  versus  Kp 
plus  the  measured  Reynolds  shear  stress  -puv  obtained  by  integrating  the  time 
series  signals  for  u(t)»v(t).  While  there  is  some  scatter,  it  may  be  seen  that 
the  variation  of  -r(y+)  is  large  enough  to  impose  an  effect  of  K  on  the  behavior 

of  uv/u2  in  the  viscous  layer.  Since  the  estimated  experimental  uncertainties 

T  _  + 
in  K  are  large,  the  uncertainties  in  the  locii  of  the  curves  T(y  )/t  are  large 
P  w 

as  well. 
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In  the  viscous  layer  the  data  increase  from  zero  towards  these  "asymptotes". 


While  the  data  for  the  lower  values  of  K  do  not  differ  from  one  another  sub- 

P 

stantially,  some  trends  may  be  discerned.  The  lowest  value  of  K  generally 

P 

leads  to  the  highest  values  of  the  Reynolds  stress  (except  for  the  first  points 


near  the  wall),  corresponding  to  the  higher  asymptotic  relation.  On  the  other 
hand,  the  higher  yields  definitely  lower  non-dimensional  Reynolds  stress  in 
the  viscous  layer,  consistent  with  its  approximate  asymptotic  curve.  (For  this 
latter  case,  is  so  large  that  -r(y+)  would  become  zero  at  y+  near  50  if  the 
convective  terms  were  indeed  negligible.) 


8.4  Skewness  and  flatness 


The  skewness  and  flatness  factors  were  calculated  from  the  time  series  for  u 
and  v  and  are  presented  as  Figures  8.5  through  8.8.  For  u  these  higher  order 
moments  are  reasonably  well  behaved  whereas  for  v  they  appear  to  suffer  a  bit 
more  experimental  scatter. 


These  factors  for  u  agree  with  the  data  of  Kreplin  [Eckelmann,  1984]  and  of 
Randolph  [1983]  for  fully-developed  flow  in  the  same  channel.  No  significant 
dependence  on  is  evident. 


The  skewness  and  flatness  data  for  t  or  (9u/3y)  show  some  variation. 

w  w 

However,  a  contributor  to  the  variation  may  be  experimental  scatter  arising 
from  low  sensitivity  of  the  sensor  in  combination  with  digital  round-off 
error.  With  the  lowest  value  of  K^,  realizations  with  the  X-probe  at  y+  ^ 
and  15  gave  S  =■  0.55  and  0.44  and  F^  a  3.16  and  2.95,  respectively. 

These  values  are  lower  than  those  of  Kreplin  (0.75  and  3.7)  which  are,  in 
turn,  lower  than  the  values  for  Su  and  F^  as  y+  approaches  zero.  At  the 
highest  pressure  gradient  K^,  two  realizations  yielded  =  0.22  and  0.24 
and  Fu  ~  2.62  and  2.62,  lower  yet. 


10 


For  the  normal  fluctuations  plotted  one  might  envision  some  trends  with 

within  the  experimental  scatter.  For  example,  Sy(y+)  may  be  shifted  a  bit  to 

the  right  for  the  highest  K  relative  to  the  lowest  in  the  viscous  layer. 

P  + 

There  is  a  general  pattern  evident  for  Sy  with  a  maximum  near  y  =10  and  a 
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zero  crossing  a  lower  y+.  Maximum  magnitudes  are  somewhat  smaller  than  for  S^, 
implying  more  symmetric  probability  density  distributions.  The  negative  values 
at  low  y+  indicate  that  large  negative  values  of  v  (sweeps?)  are  more  frequent 
than  large  positive  values  (ejections?)  near  the  wall.  There  does  not  seem  to 
be  a  clear  trend  evident  in  Fy.  However,  its  values  are  much  greater  than  F^; 
this  would  imply  that  v  frequently  takes  on  values  much  greater  than  the  mean 
(which  is  near  zero). 

8.5  Correlations 

Since  the  mean  structure  showed  only  a  gradual  variation  with  K  ,  when  it 

P 

showed  any  significant  effect  at  all,  the  rest  of  the  present  study  concen¬ 
trates  on  the  two  extremes.  Results  for  the  highest  streamwise  pressure  gradi¬ 
ent,  K  =  -0.02,  are  compared  to  those  for  K  =  -0.008,  the  lowest.  The 
P  P 

latter  is  the  condition  at  which  much  of  the  earlier  structure  results  have 

been  presented;  fully  developed  flow  at  the  maximum  flow  rate  of  the  oil 
channel  [Eckelmann,  1970,  1974;  Wallace,  Eckelmann  and  Brodkey,  1972; 

Randolph,  1983,  1988].  The  location  of  y+  =  15  was  chosen  since  it 
corresponds  to  the  most  popular  (and  most  energetic)  location  for  condi¬ 
tional  analyses  [Blackwelder  and  Kaplan,  1976;  etc.]. 

The  cross  correlation  coefficients  are  defined  as 

R 

where  u^  or  u^.  represents  any  of  the  deduced  signals  and  y+  is  the  location 
of  the  X-probe  if  involved.  Typically,  u,,  was  taken  as  the  streamwise 
velocity  fluctuation,  u.  With  i  =  j  the  result  reduces  to  the  autocorrela¬ 
tion  for  the  quantity. 

The  data  reduction  program  to  calculate  the  correlations  was  based  on  one 
written  by  Randolph  [1988].  It  was  modified  (slightly)  to  handle  four  input 
signals  and  calculate  the  "instantaneous  Reynolds  stress"  or  uv  product. 

Input  control  parameters  allowed  correlating  any  two  resulting  time  series 


♦  i  fu i<y-t> 
u.u.ly  =  i  J 

l  3  o 


u  ,  ( y , t  -  t) 


u .  u . 

l ,  rms  j  ,  rms 


dt 
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against  a  third.  All  correlations  are  presented  for  401  sample  points  at  50 

samples  per  second,  giving  a  total  period  of  about  eight  seconds  in  real  time. 

Another  short  program  was  written  to  convert  the  independent  variable  from 

■h  2 

physical  time  to  a  non-dimensional  lead  time  t  (=  tu^/v)  in  terms  of 
wall  coordinates.  While  t  has  also  been  used  as  a  symbol  to  represent 
shear  stress,  there  should  be  no  confusion  when  taken  in  context.)  Thus  the 
results  are  plotted  and  presented  as  correlation  coefficients,  non-dimension- 
alized  by  the  root  mean  square  values  of  the  fluctuations  about  their  mean 
quantities,  versus  t+  in  Figures  8.9  to  8.13. 

The  correlation  coefficients  for  fully-developed  turbulent  flow  at  y+  are 
compared  on  Figure  8.9.  The  autocorrelation  of  the  streamwise  component  peaks 

at  about  one  (=  0.99988)  as  it  should  and  shows  the  expected  symmetry.  There  does 

seem  to  be  a  slight  change  in  the  trend  of  the  slope  in  the  range  50  <  t+  <  70 

indicating  a  (very)  slightly  increased  correlation  with  some  phenomenon  at 

y+  near  60.  The  normal  component  v  shows  less  correlation  with  u,  both  in  peak 

value  and  period  of  significant  correlation.  The  u.v  correlation  is  negative, 

indicating  that  v  tends  to  be  negative  when  u  is  positive  and  vice  versa  —  a 

sweep  towards  the  wall  corresponds  to  increase  in  u  and  an  ejection  to  a 

retardation  of  u.  These  quantities  are  essentially  uncorrelated  until  t+  >  -15, 

but  the  show  some  small  level  of  correlation  to  t+  =  +  70.  The  greatest  u.v 

correlation  is  slightly  after  the  peak  in  the  u-auto  correlation,  say,  at 
+ 

t  ^  3  or  so. 

The  best  correlation  of  t  with  u  is  over  0.8  and  occurs  at  t+  =  -10.  The 

w  + 

u.t^  correlation  shows  a  greater  correlation  between  u  at  u  ==15  and  than 
between  u  at  y+  =  15  and  v  at  the  same  location!  The  peak  magnitude  is  almost 
twice  as  large  and  the  time  period  is  considerably  longer.  In  fact,  for  most 
of  the  correlation  period  the  correlation  of  u  with  is  as  large  or  larger 
than  u  with  itself  (if  the  difference  in  the  times  of  the  peak  correlation  is 
taken  into  account).  These  observations  help  explain  the  success  Chambers, 

Murphy  and  McEligot  (1983)  had  in  using  the  temporal  wall  shear  stress  in  a 
VITA  sampling  tecnnique  to  deduce  average  bursting  periods  for  laterally- 
converging  flows.  This  result  also  implies  that  it  should  be  possible  to  use 
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the  wall  sensor  as  a  detector  to  relate  the  signals  from  the  single  X-probe 
for  different  y+  locations. 

Figure  8.10  presents  the  same  correlations  for  K  =  -0.02.  Qualitatively, 

P 

the  observations  are  the  same.  Since  the  same  number  of  points  (same  real 
time)  was  used  in  the  correlations,  the  lower  value  of  u  leads  to  a  smaller 
range  of  t+.  The  individual  correlations  will  now  be  compared  directly  for 
the  two  flow  conditions  in  Figures  8.11  to  8.13;  here  the  circles  represent  the 
flow  with  the  larger  pressure  gradient  (K  =  -0.02)  and  the  square  symbols 

p 

give  the  fully-developed  results  for  reference. 


The  u-autocorrelations  are  displayed  by  Figure  8.11.  By  definition  both  peak 

at  unity  at  zero  lead,  t  +  =  0.  The  autocorrelation  at  K  =  -0.02  is  broader, 

P 

i.e.,  greater  correlation  over  a  longer  non-dimensional  time.  This  corresponds 
to  the  observation  of  Chambers,  Murphy  and  McEligot  that,  as  the  magnitude  of 
Kp  increases,  the  bursting  frequency  decreases  or  bursting  period  increases  in 
terms  of  wall  coordinates.  The  longer  average  time  between  randomly  occurring 
bursts  would  allow  eddies  to  remain  coherent  longer  in  the  vicinity  of  the 
measuring  instrument. 


Figure  8.12  compares  the.u.v  correlations.  The  peak  value  is  almost  the  same 
for  both  as  is  its  lead;  there  are  miniscule  differences  but  they  are  likely 
less  than  the  typical  experimental  uncertainty.  Again  the  correlation  is 
broader  for  the  larger  pressure  gradient.  The  fully-developed  results  are 
nearly  uncorrelated  for  t+  <  -15  whereas  the  values  for  =■  -0.02  are  just 
approaching  zero  at  the  limit  of  the  calculation  at  t+  =  -45.  After  the  peak 
u  and  v  remain  slightly  more  correlated  for  K  =  -0.02  than  for  the  fully- 
developed  run. 


The  values  of  R  at  t  =  0  reveal  consistencies  in  the  data  and  implica- 
uv 

tions  for  v'  (and  uv)  that  are  not  obvious  from  direct  examination  of  v*/u 

T 

alone.  The  quantity  R  (y  ,t  =  0)  is  sometimes  called  the  correlation  coef- 

uv 

ficient  itself;  it  may  be  calculated  directly  from  the  mean  statistics  as 
uv/u'v'  without  forming  the  more  complete  correlation.  From  Figure  8.12  we 
see  it  is  almost  exactly  the  same  for  the  two  extreme  cases  plotted  (the 
tabulated  difference  is  less  than  one  percent).  The  Reynolds  shear  stress 


-40- 


uv  decreases  about  fifty  percent  from  the  fully-developed  flow  to  the  acceler¬ 
ated  one.  This  decrease  is  shared  somewhat  equally  between  the  rms  values,  u' 
and  v',  in  yielding  the  same  value  of  R  .  However,  treating  these  results 
in  wall  coordinates  gives  a  different  picture.  The  quantity  uv/u*  decreases 
(as  expected  since  the  total  t/t^  decreases  at  this  location)  but  only  about 
14  percent.  The  corresponding  change  in  non-dimensional  rms  values  is  almost 
all  in  the  normal  component;  v'/u^  decreases  about  14  percent  and  u’/u  de¬ 
creases  less  than  one  percent.  In  other  words,  if  u’(y  J/u^  does  not  have  a 
significant  effect  of  K  ,  v’(y  )/u  must  in  order  to  remain  consistent  with 

p  T 

+  r 

the  effect  on  t(y  )/t  ;  however,  the  small  value  of  v'(<<u’)  makes  it  diffi- 

w 


cult  to  separate  this  trend  from  the  experimental  uncertainty  in  the  graph 
of  v'(y+)/uT- 


The  peak  value  of  the  correlation  between  u  and  is  approximately  the 
same  for  both  flows  (Figure  8.13)  but  it  is  slightly  later  (more  negative  lead 
time)  in  the  accelerated  case.  This  observation  may  be  a  first  indication  that 
the  inclination  of  the  fronts  of  the  sweep  events  [Kreplin  and  Eckelmann,  1979] 
changes  with  variation  of  the  pressure  gradient.  As  with  the  other  two  correla¬ 
tions  discussed,  in  this  flow  the  u.t  correlation  is  broader. 

w 


The  general  effects  of  increasing  appear  to  be  increases  in  the  correla¬ 
tions  over  longer  non-dimensional  times.  This  seems  reasonably  consistent 
with  the  idea  that  increase  of  the  bursting  period  with  K  leads  to  a 
better  organized  flow  in  the  wall  region.  Longer  non-dimensional  time  scales 
can  be  expected  to  relate  to  larger  non-dimensional  length  scales,  such  as  the 
increase  in  "linear"  layer  thickness  or  the  van  Driest  "constant"  A+  required 
in  mixing  length  predictions  of  the  mean  velocity  profile.  That  is,  the 
average  structures  must  become  larger  in  terms  of  wall  variables. 


8.6  Relation  to  measurements  with  drag-reducing  additives 


The  effects  of  polymer  injection  have  been  studied  for  several  decades 
[Hoyt  and  Fabula,  1964;  Wells  and  Spangler,  1967;  White  and  McEligot,  1970; 
Donohue,  Tiederman  and  Reischman,  1972;  Tiederman,  Luchik  and  Bogard,  1985; 
Hendr'cks,  1988,  etc.]  and  a  number  of  studies  have  examined  the  turbulence 
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structure  via  velocity  measurements,  flow  visualisation  and  wall  sensors. 

It  has  been  shown  that  polymer  additives  in  the  viscous  layer  (say  10<y+<  100) 
modify  the  structure  and  bursting  rate  (e.g.,  Tiederman,  et  al.  {1985]). 

Recently  Harder  and  Tiederman  [1989]  have  measured  fully-developed  flows 
in  two-dimensional  channels  with  low  polymer  concentration  via  laser  Doppler 
velocimetry.  Drag  reduction  of  up  to  forty  percent  was  obtained.  A  two- 
component  system  was  used  to  obtain  u,  v  and  their  statistics. 

The  turbulence  structure  quantities  deduced  by  Harder  and  Tiederman 
varied  with  the  amount  of  drag  reduction  rather  than  wall  strain  rate, 
polymer  concentration  and  channel  height,  per  se.  The  ratio  of  the  average 
period  of  the  turbulence  production  cycle  for  a  drag-reduced  flow  to  the 
average  period  for  a  water  flow  at  equal  wall  shear  stress  increases  in  a 
manner  comparable  to  the  increase  in  bursting  period  observed  with  streamwise 
pressure  gradients  by  Chambers,  Murphy  and  McEligot  [1983]  and  others. 
Likewise,  the  mean  velocity  profile  u*(y+)  increased  as  the  present  data 
do  (e.g.,  our  Figure  8.1)  but  the  magnitude  of  the  increase  appears  consider¬ 
ably  greater  for  the  polymer  flow.  The  increase  is  reminiscent  of  the  lamin- 
arizing  data  of  Perkins  and  McEligot  [1975]  and  Shehata  [1984];  this  same  ob¬ 
servation  has  been  made  by  Brodkey  earlier  11975].  The  Reynolds  shear  stress 
distribution  with  polymer  additives  (their  Figure  3.21)  shows  the  same  trend 
as  our  measurements  with  streamwise  acceleration  presented  in  Figure  8.4. 

Differences  from  the  data  of  Harder  and  Tiederman  are  evident  in  the  root 
mean  square  values  of  the  velocity  fluctuations  and  in  the  correlation  coef¬ 
ficients.  The  present  measurements  revealed  no  large  variation  of  u'/ut 

or  v’/u  with  K  .  With  typical  levels  of  drag  reduction  Harder  and  Tiederman 
t  p 

have  no  big  effect  on  v'/u  ,  particularly  in  the  viscous  layer,  but  in  their 

T 

case  with  strongly  non-Newtonian  drag  reduction  there  appears  to  be  a  de¬ 
crease  beyond  their  estimated  experimental  uncertainty.  For  u'/u^  there  is 
a  substantial  increase,  peaking  at  y+  =  30  -»  40;  this  shift  outv ard  may 
correspond  to  the  thickening  of  the  viscous  layer. 
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In  our  terms  Harder  and  Tiederman  present  Ruv(y+,-r  +  =  0)  for  correlation 

coefficients  while  our  current  data  emphasize  R.  ,(y  =  15,  i  ).Thus  we  can 

compare  ^uv(y  -  15, t  =  0).  In  our  case  we  show  no  significant  effect  of 

K  (or  K  )  whereas  Harder  and  Tiederman  have  a  decrease  of  fifty  to  near  100 
P  v 

percent  with  the  addition  of  polymers.  Theirs  is  consistent  with  the  large 
increase  in  u'(y+)/uT  and  may  be  equivalent  to  a  laminarized  flow  with  passive 
fluctuations . 

8.7  Import  for  hydroacoustic  modelling 

As  noted  above,  Thomas  and  Bull  [1983]  and  Landahl's  colleagues 
[Johansson,  Her  and  Haritonidis,  1987]  demonstrated  that  the  instantaneous 
wall  pressure  and  instantaneous  wall  shear  stress  are  related  to  one  another, 
as  indicated  earlier  by  Landahl  [1967,  1975]  who  suggested  that  unsteady 
shear  stresses  should  contribute  significantly  both  to  radiation  and  to  the 
low  wave  number  component  of  the  wall  pressure  spectrum.  Related  data  have 
been  presented  by  Zakkay,  Barra  and  Hozumi  [1979],  Dinkelacker  and 
Langeheineken  [1982;  Langeheineken ,  1981]  and  Kobashi,  Komoda  and  Ichijo 
[1986].  All  have  been  with  zero  pressure  gradients  or  fully  developed  flow 
of  air.  The  most  direct  correlation  appears  to  be  between  the  fluctuations 
of  the  normal  velocity  component  and  the  wall  pressure. 

Knowledge  of  the  turbulence  structure  is  required  to  derive  models  of  wall 
pressure  spectra  as  by  Panton  and  Linebarger  [1974],  Chase  [1980]  and  Ffowcs 
Williams  [1982],  but  no  such  models  appear  to  have  been  developed  which  in¬ 
clude  modification  of  the  turbulence  structure  by  drag  reducing  phenomena. 

Only  the  study  of  Panton  and  Linebarger  has  treated  effects  of  streamwise 
pressure  gradients. 

For  hydroacoustic  models  typical  turbulence  structure  quantities  needed 
are  source  distribution  functions  for  the  mean  shear-turbulence  component  of 
the  wall  pressure,  turbulence  -  turbulence  interaction,  etc.  The  former  may 

be  approximated  as  (dU/dy)2  •  v2 .  A  key  hypothesis  of  Panton  and  Linebarger 
[1974]  is  that  the  intensity  of  the  normal  velocity  component  scales  with  the 

Reynolds  shear  stress  and  normalized  wall  distance,  v2/uv  =  fn{y+},  but  does 
not  depend  on  pressure  gradient. 
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The  idea  that  v*  (y  )/ut  is  independent  of  streamwise  pressure  gradients 
is  tested  by  our  Figure  8.3  as  well  as  being  addressed  in  the  analysis  of 
Finnicum  and  Hanratty  [1988,  Figure  7]  and  the  direct  simulations  of  Spalart 
[1986,  19903.  A  trend  is  seen  towards  a  reduction  in  v'(y+)/u  with  increase 

T 

in  the  magnitude  of  K  ,  more  favorable  pressure  gradients.  The  question  is 

P 

how  large  is  the  effect.  For  the  range  of  the  present  data,  (-0.008<K  <-0.02, 

-6  +  p 
0<K  <2.4  x  10  )  the  reduction  is  about  thirty  percent  at  y  =25,  but  our 

comments  on  experimental  uncertainty  of  v  must  be  reiterated.  Finnicum  and 

—  6 

Hanratty  predict  a  drop  of  about  forty  percent  from  =  0  to  2.0  x  10 
while  Spalart' s  results  suggest  about  twenty  five  percent  for  increasing 
from  1.5  x  10  6  to  2.5  x  10  6,  both  at  y+  =  25. 


9.  CONCLUDING  REMARKS  AND  RECOMMENDATIONS 


To  provide  bases  for  incorporating  the  effects  of  streamwise  pressure 
gradients  into  hydroacoustic  models  for  the  wave  vector-frequency  spectrum  of 
wall  pressure  fluctuations  in  turbulent  boundary  layers,  the  turbulence 
structure  has  been  studied.  Efforts  were  concentrated  on  the  viscous  layer 
which  often  dominates  momentum,  energy  and  mass  transfer  from  the  wall.  (Here 
we  define  the  viscous  layer  with  Bradshaw  [1969]  as  the  region  of  a  turbulent 
flow  where  viscous  effects  are  significant  though  not  necessarily  dominant, 
typically  0  <  y+  <  30  in  a  fully-developed  flow.) 

Earlier  measurements  of  the  temporal  wall  shear  stress  by  Chambers,  Murphy 
and  McEligot  [1983]  in  a  laterally-converging  flow  had  led  to  the  conclusions 
that 

o  non-dimensional  bursting  frequencies  decrease  greatly  as  the  acceleration 
parameter  increases 

o  the  time  history  of  the  conditionally-averaged  wall  shear  stress  was  less 
affected 

o  results  presented  in  terms  of  wall  scaling  were  less  sensitive  to  accelera¬ 
tion  than  those  normalized  by  outer  scaling. 

These  observations  came  from  experiments  in  an  air  flow  where  the  half-inch 
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(about  13mm)  plate  spacing  and  thin  shear  layer  prohibited  fluid  velocity 
measurements . 

The  present  studies  were  conducted  in  the  22cm  wide  Reichardt/Eckelmann  oil 
channel  at  the  Max  Planck  Institut  fur  Strbmungsf orschung  to  obtain  good 
spatial  and  temporal  resolution.  Simultaneous  time  series  data  were  obtained 

with  an  X-probe  and  a  wall  sensor  to  determine  u(t),  v(t)  and  t  (t).  Measure- 

w 

ments  with  the  X-probe  at  y  =  5,  7,  10,  15,  25  and  the  centerplane  concen¬ 
trated  on  four  sets  of  conditions: 


0 

K  - 

-0.02, 

K 

-  2.4 

x  10  6  , 

Re 

=  5600 

P 

V 

s 

0 

K  = 

-0.011, 

K 

s  1.6 

x  10  6  , 

Re 

=  8300 

P 

V 

s 

o 

K  - 

-0.011, 

K 

-  o, 

Re 

=  5600 

P 

V 

s 

0 

K  - 

-0.008, 

K 

*  0, 

Re 

=  7400 

P 

V 

s 

In  terms  of  wall  layer  scaling,  the  following  effects  were  found  when  the 
(non-dimensional)  streamwise  pressure  gradient  increased: 

o  the  mean  velocity  profile,  u+(y+),  increased,  corresponding  to  a 

thickening  of  the  layer  dominated  by  viscous  effects 

o  at  y+  -  15,  correlation  coefficients,  R  (t  +  ),  R  (t  +  )  and 

uu  uv 

R  (t  )  broadened,  consistent  with  the  decrease  in  bursting 

UT 

frequencies  observed  by  Chambers,  Murphy  and  McEligot 

The  following  aspects  of  the  turbulence  structure  showed  no  significant 
effects  in  the  viscous  layer: 

o  u'(y+)/u^,  root  mean  square  of  the  streamwise  velocity  fluctuations 

o  Su(y+),  skewness  factors  of  the  streamwise  velocity  fluctuations 

o  F  (y  ),  flatness  factors  of  the  streamwise  velocity  fluctuations 
u 

o  maximum  amplitudes  of  the  correlation  coefficients  R  (t  )  and  R  (t  )  (and 

uv  UT 

R  by  definition)  at  y  ^15 

Experimental  uncertainties  in  determining  the  instantaneous  normal  velocity  com¬ 
ponent  v  (i.e.,  calculation  of  a  small  difference  between  larger  quantities) 
hamper  conclusions  on  the  statistics  of  v.  But  it  appears  a  tentative  con- 
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elusion  is  warranted  that  there  is  no  major  effect  on  S  (y+)  or  F  (y+). 

v  v 

Calculations  of  the  non-dimensional  spectra  are  in  process.  It  is  recom¬ 
mended  that  they  be  completed.  Modification  of  the  versatile  pattern  recog¬ 
nition  program  of  Brodkey,  to  correspond  to  the  details  of  the  present  measure¬ 
ments,  has  also  been  initiated.  In  addition  to  pattern  recognition,  this 
program  applies  the  VITA  technique  [Blackwelder  and  Kaplan,  1976]  with  a 
variety  of  detection  schemes  suggested  in  the  literature  plus  a  number  of 
quadrant-splitting  approaches  [Wallace,  Eckelmann  and  Brodkey,  1972],  It  is 
likewise  recommended  that  this  aspect  of  the  study  be  completed  to  determine 
further  in  which  ways  streamwise  pressure  gradients  affect  the  turbulence 
structure  in  the  dominant  viscous  layer. 

These  data  and  the  comparable  measurements  of  Randolph  [1983]  show  signs  of 
interaction  between  the  probe  and  the  wall  as  the  probe  approaches  the  wall.  It 
is  believed  that  relative  or  comparative  results  are  not  significantly  affected 
although  the  mean  wall  shear  stress  and  the  mean  flow  angle  vary  somewhat. 
Therefore,  it  is  recommended  that  these  data  be  supplemented  by  non-intrusive 
measurements  with  laser  Doppler  velocimeters  in  the  oil  channel.  Such  instru¬ 
mentation  could  offer  the  further  benefit  of  more  direct  measurement  of  the 
normal  velocity  and,  possibly,  more  accurate  determination  of  its  statistics. 
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2  Conditionally-averaged  wall  shear-stress  fluctuations  in  laterally 
converging  flow;  11,000  <  Re^  <  13,800;  0  <  Ky  <  3.4  x  10 
[Chambers,  Murphy  and  McEligot,  1983]. 
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3  Average  bursting  periods  (inverse  of  bursting  frequency)  for  laterally 
converging  flows  [Chambers,  Murphy  and  McEligot,  1983]. 
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Fig.  7.1  Mean  velocity  predictions  using  function  A(K  )  deduced  from  Huffman  and 
Bradshaw  11972]  in  van  Driest  mixing  length  model  [1956). 
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Fig.  7.2  Influences  of  measurements  of  X-probe  location  and  estimates  of  pressure 
gradient  parameter  [McEligot,  1984]. 
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Fig.  7.4  Graphical  representation  of  correction  process  for  effect  of 
fluid  temperature. 
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Fig. 


Root  mean  square  values  of  streamwise  velocity  component  fluctuations 
Symbols  as  in  Pig.  8.1. 
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Fig.  8.3  Root  mean  square  values  of  normal  velocity  component 
fluctuations.  Symbols  as  in  Fig.  8.1. 


30 


Fig.  8.4  Variation  of  Reynolds  shear  sv  ss. 
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Distribution  of  skewness  factors  of  streamwise  velocity  fluctuations 
Symbols  as  in  Fig.  8.1. 


Fig.  8.6  Distribution  of  skewness  factors  of  normal  velocity  fluctuations 
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Distribution  of  flatness  factors  of  normal  velocity  fluctuations. 


Fully  Developed 
Re. =  7400 


correlation  coefficients  for  fully-developed  flow,  K  =  -0.008,  y  =  15. 


Comparison  of  autocorrelation  coefficients  of  streamwise  fluctuations. 
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Appendix  A.  MEAN  STATISTICS  OF  VISCOUS  LAYER 

’•Fully-developed"  flow 


K  =  -0.00831,  u 

P  t 

=1.05  cm/s,  pump 

speed  = 

1790  rpm 

Approximate  y+ 

7 

10 

15 

25 

Centerplane 

T,  °C 

25.02 

25.0 

24.98 

24.95 

24.72 

v,  cm2/sec 

0.05977 

0.05980 

0.05984 

0.05988 

0.06024 

y  ,mm 

3.224 

4.922 

7.529 

12.702 

109.3 

u,  cm/sec 

5.79 

8.12 

10.91 

14.04 

20.32 

v,  cm/sec 

-0.295 

-0.314 

-0.200 

-0.0499 

0.0861 

u ' ,  cm/sec 

2.06 

2.71 

3.04 

2.79 

0.884 

v' ,  cm/sec 

0.222 

0.290 

0.408 

0.615 

0.622 

uv  (cm/sec)2 

-0.244 

-0.352 

-0.450 

-0.607 

-0.0404 

su 

0.606 

0.331 

-0.0673 

-0.432 

-0.614 

sv 

0.179 

0.316 

0.286 

-0.0228 

0.139 

FU 

2.97 

2.43 

2.16 

2.52 

3.63 

FV 

4.34 

5.31 

6.28 

4.64 

3.93 

y+ 

5.64 

8.60 

13.15 

22.17 

189.6 

u+ 

5.54 

7.77 

10.43 

13.43 

19.44 

v+ 

-0.282 

-0.301 

-0.191 

-0.0477 

0.0824 

u’  + 

1.98 

2.59 

2.91 

2.67 

0.846 

v'  + 

0.214 

0.278 

0.390 

0.589 

0.595 

uv/u  2 

T 

-0.223 

-0.322 

-0.412 

-0.555 

-0.0370 

v/u 

-0.0509 

-0.0387 

-0.0183 

-0.00356 

+0.00424 

u '  /u 

0.356 

0.334 

0.279 

0.199 

0.0435 

< 

' — 

C  1 

0.0384 

0.0358 

0.0374 

0.0438 

0.0306 
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Appendix  A.  MEAN  STATISTICS  OF  VISCOUS  LAYER  (Cont’d. 

"Fully-developed"  flow 


K  =  -0.0111,  u 
P  t 


0.801  cm/sec,  pump  speed  =  1350  rpm 


Approximate  y+ 

5 

7 

10 

15 

25 

T,  °C 

24.9 

24.9 

24.82 

24.82 

24.85 

v,  cm2/sec 

0.05996 

0.05996 

0.06009 

0.06009 

0.06004 

y  ,mm 

3.831 

5.585 

8.136 

12.001 

19.717 

u,  cm/sec 

4.02 

5.54 

7.49 

9.41 

11.41 

v,  cm/sec 

-0.203 

-0.224 

-0.195 

-0.123 

-0.0195 

u’ ,  cm/sec 

1.47 

1.89 

2.19 

2.30 

1.96 

v' ,  cm/sec 

0.157 

0.179 

0.235 

0.325 

0.483 

uv  (cm/sec)2 

-0.165 

-0.190 

-0.214 

-0.285 

-0.345 

su 

0.653 

0.398 

0.0988 

-0.210 

-0.521 

sv 

-0.105 

0.305 

0.409 

0.0501 

-0.0518 

Fu 

3.13 

2.55 

2.22 

2.21 

2.66 

Fv 

3.66 

4.68 

6.39 

5.95 

4.18 

y+ 

5.12 

7.46 

10.85 

16.00 

26.31 

u+ 

5.02 

6.92 

9.36 

11.74 

14.25 

v+ 

-0.253 

-0.280 

-0.244 

-0.154 

-0.0243 

u'  + 

1.84 

2.36 

2.74 

2.87 

2.45 

v*  + 

0.196 

0.224 

0.294 

0.406 

0.603 

uv/u2 

-0.257 

-0.296 

-0.333 

-0.445 

-0.537 

v/u 

-0.0504 

-0.0405 

-0.0261 

-0.0131 

-0.0017 

u '  /u 

0.366 

0.341 

0.293 

0.245 

0.172 

v'  /u 

0.0391 

0.0324 

0.0314 

0.0346 

0.0423 

Centerplane 

24.74 

0.06021 

109.4 

15.39 

0.0798 

0.701 

0.468 

-0.0326 

-0.688 

+0.111 

3.79 

4.04 

145.5 
19.22 

0.0996 

0.875 

0.584 

-0.0407 

+0.00518 

0.0455 

0.0304 
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Appendix  A.  MEAN  STATISTICS  OF  VISCOUS  LAYER  (Cont'd.) 


Laterally  converging  (spatially  accelerated)  flow 


K  =  -0.0108,  u 

=  1.16  cm/s,  pump  speed  -  1790 

rpm 

P  T 

Approximate  y+ 

7 

10 

15 

25 

T,  °C 

25.07 

25.05 

25.0 

24.98 

v,  cm2/sec 

0.05970 

0.05973 

0.05980 

0.05984 

y  ,mm 

3.50 

5.20 

7.801 

13.001 

u,  cm/sec 

7.46 

10.17 

13.10 

16.32 

v,  cm/ sec 

-0.379 

-0.386 

-0.289 

-0.178 

u  * ,  cm/sec 

2.51 

3.09 

3.32 

2.93 

v' ,  cm/sec 

0.341 

0.391 

0.471 

0.615 

uv  (cm/sec)2 

-0.427 

-0.486 

-0.597 

-0.666 

su 

0.544 

0.220 

-0.124 

-0.468 

Sv 

-0.0293 

0.270 

0.356 

0.0308 

fu 

2.78 

2.28 

2.17 

2.57 

Fy 

3.49 

3.96 

5.08 

4.60 

y+ 

6.80 

10.10 

15.13 

25.20 

u+ 

6.43 

8.77 

11.30 

14.07 

v+ 

-0.327 

-0.333 

-0.249 

-0.153 

u'  + 

2.16 

2.66 

2.86 

2.53 

v*  + 

0.294 

0.337 

0.406 

0.530 

uv/uf 

-0.318 

-0.361 

-0.444 

-0.495 

v/u 

-0.0508 

-0.0380 

-0.0221 

-0.0109 

u '  /u 

0.337 

0.304 

0.253 

0.180 

V'  /u 

0.0457 

0.0384 

0.0360 

0.0377 

Centerplane 

24.72 

0.06024 

110.0 

22.37 

-0.191 

0.822 

0.517 

-0.103 

-0.697 

0.150 

3.79 

4.06 

211.8 

19.29 

-0.165 

0.709 

0.446 

-0.0766 

0.00854 

0.0367 

0.0231 
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Appendix  A.  MEAN  STATISTICS  OF  VISCOUS  LAYER  (Cont'd. 


Laterally 

converRinR 

(spatially 

accelerated)  flow 

K  =  -0.0201,  u 

P  t 

=  0.793  cm/s,  pump  speed  =  1200 

rpm 

Approximate  y+ 

5 

7 

10 

15 

25 

T,  °C 

25.08 

25.0 

25.05 

25.08 

25.09 

v,  cm2/sec 

0.05968 

0.05980 

0.05973 

0.05968 

0.05966 

y  ,mm 

3.84 

5.592 

8.121 

12.019 

19.719 

u,  cm/sec 

3.  75 

5.39 

7.36 

9.58 

11.75 

v,  cm/sec 

-0.213 

-0.215 

-0.187 

-0.139 

-0.0917 

u ' ,  cm/sec 

1.21 

1.  79 

2.15 

2.29 

1.96 

v' ,  cm/sec 

0.165 

0.195 

0.219 

0.2680 

0.347 

uv  (cm/sec)2 

-0.125 

-0.186 

-0.186 

-0.225 

-0.262 

su 

0.757 

0.563 

0.205 

-0.197 

-0.551 

sv 

0.437 

-  0.0506 

0.470 

0.433 

0.0559 

Fu 

3.37 

2.69 

2.24 

2.15 

2.64 

f'v 

4.04 

3.57 

5.33 

5.09 

5.20 

y+ 

5.11 

7.43 

10.80 

15.99 

26.25 

u+ 

4.72 

6.793 

9.27 

12.06 

14.80 

v+ 

-0.268 

-0.271 

-0.235 

-0.175 

-0.116 

u,+ 

1.53 

2.26 

2.71 

2.88 

2.47 

V>  + 

0.208 

0.246 

0.276 

0.338 

0.437 

uv/u| 

-0.200 

-0.234 

-0.295 

-0.356 

-0.416 

v/u 

-0.0568 

-0.0399 

-0.0254 

-0.0145 

-0.00780 

u  *  /u 

0.324 

0.332 

0.292 

0.239 

0.167 

V’  /u 

0.0440 

0.0362 

0.0298 

0.0280 

0.0296 

Center-plane 

25.05 

0.05973 

110.0 

15.25 

-0.0280 

0.572 

0.336 

-0.0349 

-0.942 

0.0824 

4.75 

4.07 

146.3 

19.20 

-0.0352 

0.720 

0.423 

-0.0553 

-0.00184 

0.0375 

0.0220 


-74- 


Appendix  B.  MEAN  VELOCITY  PROFILES 

"Fully-developed"  flow 


K  =  -0.00831,  pump  speed  =  1790  rpm,  V  =  17.23  cm/sec 

p  4  b 

Re  =  2V.s/\>  -  1.25  x  10  ,  Re  =  u  s/v  =7410. 


b 

.  s 

c 

+ 

+ 

y  (mm) 

u ( cm/s ) 

u 

_ 1 _ 

0 

0 

0 

0 

3.773 

5  .  74 

5.49 

6.56 

6.492 

9.29 

8.88 

11 . 29 

9.556 

11.94 

11.42 

16.61 

14.556 

14.25 

13.62 

25.31 

19.489 

15.23 

14.56 

33.88 

29.508 

16.87 

16.13 

51.30 

39.3 

17.53 

16 . 76 

68.3 

49.3 

18.24 

17.43 

85.7 

59.3 

18.78 

17.96 

103.1 

69.3 

19 . 36 

18.51 

120.5 

89.3 

20.00 

19.12 

155 . 3 

109 . 3 

20.39 

19.49 

190.0 
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K 

P 

Re 


Appendix  B,  MEAN  VELOCITY  PROFILES  (Cont'd.) 

"Fullv-developed"  flow 


1 

o 

o 
►— i 
h-1 

pump  speed  =  1350  rpm. 

Vv  =  12.82 
b 

cm/sec , 

V  2s/v  = 

-  9320,  Re  =  u  s/v  = 

5600. 

b 

s  c 

+ 

y  (mm) 

u(cm/s) 

u 

.  Y _ 

0 

0 

0 

0 

3.805 

3.79 

4 . 73 

5.07 

6.477 

6.01 

7.49 

8.63 

8.964 

7.53 

9.40 

11.94 

14.717 

10.26 

12.81 

19.60 

19.241 

11.04 

13.78 

25.63 

29.099 

12.54 

15.64 

38.76 

39.4 

13.15 

16.42 

52.5 

49.4 

13.69 

17.08 

65.8 

59.4 

14.16 

17.67 

79.1 

69.4 

14.45 

18.04 

92.4 

89.4 

15.12 

18.87 

119.1 

109.4 

15.39 

19.20 

145.7 
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Appendix  B.  MEAN  VELOCITY  PROFILES  (Cont’d.) 


Laterally  converging  (spatially  accelerated)  flows 


K  -  -0.0108,  V.  =  19.48  cm/sec, 
p  b 


Re  =  V.  2s/ v  - 
D 

1.43  x  104,  Re 

s 

=  u  s / v  —  8290. 
c 

+ 

+ 

y  (mm) 

u  ( cm/s ) 

u 

_ 

0 

0 

0 

0 

5.000 

9.39 

8.09 

9.64 

10.000 

14 . 38 

12.39 

19.28 

14.020 

16 . 76 

14.44 

27.03 

27.7 

19.04 

16.40 

53.4 

55.2 

20.95 

18.04 

106.4 

82.7 

22.07 

19.01 

159.5 

110.2 

22.64 

19.50 

212.5 

Kp  *  -0.0201, 
Re  =  Vfc2s/v  = 

V]-,  =  12.88  cm/sec, 
9530,  Res  =  ucs/v  = 

5600. 

+ 

+ 

y  (mm) 

u(cm/s) 

u 

-JL _ 

0 

0 

0 

0 

5.339 

4.99 

6.29 

7.09 

10.339 

8.50 

10.71 

13.73 

14.320 

10.49 

13.21 

19.02 

28.1 

12.56 

15.81 

37.3 

55.6 

14 . 14 

17.81 

73.8 

83.1 

14.92 

18.79 

110.4 

110.6 

15.23 

19.1  P 

146.9 
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Appendix  C.  Computer  program  TAU2.FOR 


As  described  in  Section  7,  the  technique  for  determination  of  the  mean  wall 
shear  stress,  from  a  single  mean  velocity  measurement  and  the  pressure  gradient, 
was  modified  to  account  for  the  convective  terms  in  the  streamwise  momentum 
equation.  To  accomplish  the  calculations,  Fortran  program  TAUW.FOR  tMcEligot, 
1984,  Appendix  B]  was  revised  to  TAU2.F0R  which  is  provided  here.  A  listing  of 
the  program  follows,  and  then  copies  of  the  interactive  iterative  input  state¬ 
ments  and  an  example  of  the  typical  output.  The  bases  of  the  program  are 
described  in  Section  7. 


Program  listing 


VAX  FORTRAN  V4 . 8-276  Page  1 

DI SK$D1 SK2 : [ LEHMANN . OLDVAX ] TAU2 . FOR ; 9 


19-Mar-1990  18:08:31 
19-Mar-1990  17:39:33 


0001 

CC 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

C 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  D.  M.  McELIGOT  bei  H.  ECKELMANN  19  Mar  1990 

C  1745 

C  File  TAU2.FOR 

C 

C  Calculation  of  wall  shear  stress  with  acceleration  correction 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  PROGRAM  TO  DEDUCE  WALL  SHEAR  STRESS 
C  WITH  VAN  DRIEST  MODEL 

C  HUFFMAN  &  BRADSHAW  VERSION  (JFM,  1972) 

C 

COUT  DOUBLE  PRECISION  Y1,Y2,  Y3,  FI , F2 , F3 , E,G, K , YL,  YM,  U1,U2,U3,  KP,D 

IMPLICIT  INTEGER* 4  (I-N),  REAL*8(A-H,0-Z) 

REAL*8  K,  KP 
CHARACTER *9  TAG,  ZEIT*8 
DIMENSION  Cl (1000) 

C 

K  -  0.4 
YL0  -=26.0 
A  -  0.0 
IP  -  2 
ICOR  -  0 
C 

CALL  TIME(ZEIT)  !  Find  system  time 

CALL  DATE ( TAG )  !  Find  system  date 

WRITE( 7 , * ) 

WRITE( 7 ,11102 )  ZEIT,  TAG 

11102  FORMAT!/  '  Start  of  calculation  at 

+  A8, '  on  ' ,A9,/  ) 

WRITE(7, 11101) 

11101  FORMAT!/ 

+  »  *********************************************************** < ^ / 

+  '  D.  M.  McELIGOT  bei  H.  ECKELMANN  '// 

V  ********************************************************** v/) 

C 

10  WRITE! 7,*)  'DATA  REDUCTION  FOR  WALL  SHEAR  STRESS' 
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Program  listing  (Conl'd.) 


0040 

WRITE! 6,*)  'TEMP(C)?  TOLERANCE?' 

0041 

READ)  5,  *  )  T,  EPS 

0042 

WRITE (6,*)  'DZ  PROBABLY  NEGATIVE' 

0043 

WRITE! 6 , * )  ' DZ (MM  OIL)?  DX(CM)?' 

0044 

READ! 5, * )  DZ,  DX 

0045 

C 

0046 

C 

OIL  PROPERTIES  (REICHARDT, 1953 ;ECKELMANN, 1970 ) 

0047 

RO  =  0.8358  -  0.000685*(T  -  20.) 

0048 

C 

GM/CC 

0049 

XMU  -  4.99  -  0.13375MT  -  24.91) 

0050 

C 

CENTI POISE 

0051 

XMU  -  0 . 001 *XMU 

0052 

C 

KG/M. SEC 

0053 

XNU  -  10 . *XMU/RO 

0054 

C 

CM2/SEC 

0055 

c 

0056 

DZ  *  0.1 *DZ 

0057 

c 

CM  OIL 

0056 

WRITE(7,*)  '  DZ  -  ' , DZ , '  CM,  DX  -  ',DX,'  CM' 

0059 

WRITE! 7,*)  '  TEMP  -  ',T,'  C,  RHO  -  ',RO,'  GM/CM3 ' 

0060 

WRITE! 7 , * )  '  VISC  -  ' ,XMU, '  KG/M. SEC,  KIN  VISC  - 

0061 

+  '  CM2/SEC' 

0062 

c 

0063 

c 

START  LOOP  FOR  NEW  DATA  POINT 

0064 

c 

0065 

30 

WRITE! 6,*)  '  Apply  acceleration  correction  (Y-l.N-1 

0066 

READ! 5,*)  ICOR 

WRITE!  6,*)  '  YKMM)  ?  UI(  CM/SEC)  ?  N  ? ' 

0067 

0068 

READ! 5,*)  YI,  UI,  N 

0069 

WRITE! 6 , * )  '  Estimate  of  u,cl(cm/sec)  ?  >' 

0070 

READ! 5,*)  UCL 

0071 

WRITE! 7 , * )  ’  ICOR-  ' , I COR , ' ,  u , cl ( cro/sec )  -  ' ,UCL 

0072 

J  -  1 

0073 

YI  -  0.1*YI 

0074 

C 

CM 

0075 

DCM  -  (YI  -  A)/N 

0076 

UICM  -  0.0 

0077 

UT2M  -0.0 

0078 

IPRO  -  0 

0079 

XCR  -  ICOR 

0080 

DO  35  I  -  1,  N 

0081 

CHI)  -  0.0D00 

0082 

35 

CONTINUE 

0083 

C 

0084 

C 

FIRST  ESTIMATE 

0085 

C 

0086 

UT2  -  1.  -  (0. 5*980. 7*DZ*YI*YI )/(XNU*DX*UI ) 

0087 

UT2  -  UT2*XNU*UI/YI 

,  XNU, 


)?  >' 
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0088 


C 


0089 

C 

START  OF  ITERATION  LOOP 

0090 

c 

0091 

40 

CONTINUE 

0092 

UT  -  SQRT(UT2) 

0093 

UT3  -  UT*UT2 

0094 

0095 

0096 

0097 

KP  -  980 . 7*XNU*DZ/( UT3*DX ) 

IFIKP.LT. (-0.0249) )  KP  -  -0.0240 

CF  m  2.0  *  UT2  /  UCL**2 

D  -  DCM*UT/XNU 

0098 

YM  -  YI*UT/XNU 

0099 

IC  -  0 

0100 

A  -  0.0 

0101 

o 

o 

M 

H 

0102 

Ul  -  0.0 

0103 

FI  -  1.0 

0104 

c 

0105 

c 

APPROXIMATION  OF  HUFFMAN  A+ 

0106 

YL  -  YL0  -  31 . 3*LOG( 1 .  +  (KP/0.025)  ) 

0107 

c 

0108 

UO  -  0.0 

0109 

U02  -  0.0 

0110 

WRITE( 7 , 20 )  K,  YL,  KP,  YM,  N,  IP,  D 

0111 

20 

FORMAT ( IX , //I OX , ' MEAN  VELOCITY  PROFILE 

0112 

1  10X, 'VAN  DRIEST  (HUFFMAN/BRADSHAW)  1 

0113 

2  5X, 'KAPPA  -  ' , F9.6 , '  Y+,L  -  ',F9.4 

0114 

3  5X,'Y+,MAX  -  ' ,0PF9.4, '  N  -  ',13,' 

0115 

4  '  DY+  •  ' , 1PE12 . 4,/// 

0116 

5  6X, 'Y+' ,7X, 'U+' ,6X, ' (U+  -  Y+)/U+'r 

0117 

6  2X,'(U+  -  Y+)/Y+' ,2X, ' (U+  -  UC)/U+' 

0118 

WRITE( 7,22)  Yl,  Ul ,  UO,  U02 

0119 

N1  -  N  +  1 

0120 

NMl  -  N  -  1 

0121 

DO  55  I  -  1,NM1 ,2 

0122 

1C  -  IC  +  1 

PRINT  INT 


1PE12.4/ 
-  ',  13, 
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Program 

listing. 

(Cont’d. ) 

0123 

Y2  -  Y1  +  D 

0124 

Y3  -  Y2  +  D 

0125 

C 

0126 

C 

VAN  DRIEST  MODEL  USED 

0127 

E  -  K*Y2* ( 1 .  -  EXP(-Y2/YL) ) 

0128 

E  -  E*E 

0129 

0130 

0131 

0132 

+ 

C234567 

G  «  1.  4  RP*Y2 

* ( 1 . 0- ( XCR*CF*CI (I+1)/(2.0*Y2)  )  ) 

F2  -  (-1.  4  SORT ( 1 .  4  4 . *E*G ) )/( 2 . *E ) 

0133 

E  -  K*Y3* ( 1 .  -  EXP(-Y3/YL)) 

0134 

E  -  E*E 

0135 

0136 

0137 

4 

G  -  1.  +  KP*Y3 

* ( 1 . 0-( XCR*CF*CI (I+2)/(2.0*Y3)  )  ) 

F3  -  (-1.  4  SQRT(1.  +  4 . *E*G) )/( 2 . *E ) 

0138 

C 

0139 

c 

SIMPSON'S  RULE  INTEGRATION 

0140 

c 

0141 

U3  -  ( D* ( FI  4  4 . *F2  4  T3)/3.)  4  ul 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

c 

c 

c 

c 

50 

CORRECTION  FOR  CONVECTIVE  TERMS  IN  ACCELERATED  FLCWS 

IF(ICOR.NE.l)  GO  TO  50 

U2  -  (Ul  4  U3 )  /  2.0 

CI(l4l)  -  Cl ( I )  4  (Ul*Ul4U2*U2)*D/2.0 

CI(I42)  -  CI(l4l)  4  (U2*U24U3*U3)*D/2.0 

CONTINUE 

Ul  -  U3 

0152 

Yl  -  Y3 

0153 

Fl  «  F3 

0154 

IF( IPRO.NE. 1 )  GO  TO  55 

0155 

IF( IC.NE. IP)  GO  TO  55 

0156 

IC  -  0 

0157 

UO  -  (U3  -  Y3)/U3 

0158 

U02  -  (U3  -  Y3 )/Y3 

0159 

UC  -  Y3*(l.  4  0 . 5*KP*Y3 ) 

0160 

U03  -  (U3  -  UC)/U3 

0161 

WRITE( 7 , 22 )  Y3,  U3,  UO,  U02,  U03 

0162 

22 

FORMAT ( 2X,F8 . 3, 2X,F8 . 3, 2X, 1PE12 . 4 , 3( 2X,E12 . 4 )  ) 
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Program  listing  (Cont’d.) 


0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 

0181 

0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 


55  CONTINUE 

UO  -  (U3  -  Y3)/U3 

U02  -  (U3  -  Y3)/Y3 

UIC  -  U3*UT 

ERR  «  (UIC  -  UIJ/UI 

WRITE( 7,22)  Y3,  U3,  UO,  U02 

WRITE(6,*)  YI ,UI ,UT,Y3,U3,UIC 

WRITE! 6 , * )  '  Convergence  error ( fraction)  -  ' , ERR 
WRITE(6,*)  '  Trial  K,p  -  KP 

WRITE ( 7 , * )  '  ITERATION:  y,meas( cm) ,  u ,meas ( cm/sec ) ,  u, tau( cm/sec > ' 
WRITE! 7,*)  YI,UI,UT 

WRITE! 7 , * )  '  y,+(pred) ,  u+(pred),  UIC' 

WRITE! 7 , * )  Y3,U3,UIC 

WRITE! 7 , * )  '  Convergence  error! fraction)  -  ' ,ERR 

WRITE! 6,*)  'MORE?  (1-YES)  '  !  Iteration  control 

READ! 5,*)  IA 

UT2P  -  ( (UI-UIC)*(UT2-UT2M)/(UIC-UICM) )  +  UT2 
UICM  -  UIC 
UT2M  -  UT2 
UT2  -  UT2P 

IF(IA.EQ.l)  GO  TO  40 
UT  -  SQRT(UT2) 

UT3  -  UT*UT2 

CF  -  2.0  *  UT2  /  UCL**2 

KP  -  980 . 7*XNU*D7./(  UT3*DX ) 

WRITE! 7, 11103) 

11103  FORMAT!/'  •****************•*******************************'/) 

WRITE! 7 , * )  '  ICOR  -  ' , I COR , ' ,  u, tau( cm/sec )  -  ',UT 

WRITE! 7 , * )  '  C,f(based  on  u,cl)  -  ',CF, 

+  K, p  -  KP 

WRITE! 7,*)  '  (ICOR  -  1  -  acceleration  correction  included' 
WRITE! 7,*)  ’  -  0  -  not  included  )' 

WRITE(7, 11103) 

WRITE! 6,*)  'ANOTHER  POINT?  < 1-YES) ' 

READ! 5 , * )  IB 

IF( IB.EQ.l )  GO  TO  30 

WRITE! 6,*)  'NEW  T  OR  DP  ?  (1-YES)' 

READ! 5, * )  ID 

IF(ID.EQ.l)  GO  TO  10 
WRITE! 7, 11102)  ZEIT,  TAG 
WRITE! 7, 11101) 

C 

STOP 

END 
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Typical  interactive  input  from  terminal 


Unit  6 

provides  the 

prompts  to  the  terminal  screen  and  Unit  5  is  the  user 

response  at  the  terminal. 

0039  10 

WRITE(7,*) 

'DATA  REDUCTION  FOR  WALD  SHEAR  STRESS' 

0040 

WRITE ( 6  ,  * ) 

'TEMP(C)?  TOLERANCE?' 

0041 

R£AD( 5 , * ) 

T,  EPS 

0042 

WRITE!  6, * ) 

'DZ  PROBABLY  NEGATIVE' 

0043 

WRITE ( 6 , * ) 

'  DZ  ( MM  OIL)?  DX ( CM ) ? ' 

0044 

R£AD( 5 , * ) 

DZ,  DX 

0063 

c 

START 

LOOP  FOR  NEW  DATA  POINT 

0064 

c 

0065 

30 

WRITE( 6 , * ) 

'  Apply  acceleration  correction  (Y-1,N«0)?  >' 

0066 

READ( 5 , * ) 

I  COR 

0067 

WRITE* 6,*) 

'  YI(KM)  ?  UI ( CM/SEC )  ?  N  ? ' 

0068 

READ(5,*) 

YI,  UI,  N 

0069 

WRITE(6,*) 

'  Estimate  of  u,cl(cm/6ec)  ?  >’ 

0070 

read; 5,*) 

UCL 

0169 

write; 6,*) 

YI ,UI,UT,Y3,U3,UIC 

0170 

write; 6,*) 

'  Convergence  error ( fraction)  -  ' ,ERR 

0171 

WRITE(6,*) 

'  Trial  K,p  -  ' ,  KP 

0177 

WRITE ( 6 , * )  'MORE?  (1-YES)  ' 

1  Iteration  control 

0178 

read; 5,*)  ia 

0196 

WRITE* 6,*)  'ANOTHER  POINT?  (1-YES)' 

0197 

read; 5, * )  ib 

0198 

IF(IB.EQ.l)  GO  TO  30 

0199 

WRITE( 6 , * )  'NEW  T  OR  DP  ?  (1-YES)' 

0200 

READ( 5 , * )  ID 

0201 

IF(ID.EQ.l)  GO  TO  10 

0202 

WRITE(7, 11102)  ZEIT,  TAG 

0203 

WRITE(7, 11101) 

0204 

C 

0205 

STOP 

0206 

END 
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Typical  output  from  interactive  session 


DATA  REDUCTION  FOR  WALL  SHEAR  STRESS 
DZ  =  -3 . 0000000000000000E-02  CM,  DX  =  174.7000000000000  CM 

TEMP  =  25.08000000000000  C,  RHO  =  0.8323202000000000  GM/CM3 

VISC  -  4.9672625000000000E-03  KG/M. SEC,  KIN  VISC  -  5 . 9679706199609238E-02 

CM2/SEC 

ICOR  =  1,  u, cl (cm/sec)  =  15.25100000000000 


MEAN  VELOCITY  PROFILE 

VAN  DRIEST  (HUFFMAN/BRADSHAW)  MODEL 

KAPPA  =  0.400000  Y+,L  =  104.1371  K,P  -  -2.2940E-02 

Y+,MAX  =  15.2958  N  =  150  PRINT  INT  -  2  DY+  =  1.0197E-01 


Y+  U+  (U+  -  Y+)/U+  (U+  -  Y+)/Y+  (U+  -  UC)/U+ 

0.0C0  0.000  0. 0000E+00  0 . OOOOE+OO 

15.296  11.833  -2.9268E-01  -2.2641E-01 

ITERATION:  y,meas(cm),  u,meas( cm/sec ) ,  u, tau( cm/sec ) 
1.201900000000000  9.579000000000000  0.7595033691333739 

y,+(pred),  u+(pred),  UIC 

15.29577066462467  11.83262820706867  8.986920988971251 

Convergence  error ( fraction)  =  -6 . 1810106590327678E-02 


MEAN  VELOCITY  PROFILE 

VAN  DRIEST  (HUFFMAN/BRADSHAW)  MODEL 

KAPPA  =  0.400000  Y+,L  -  82.1830  K,P  -  -2.0847E-02 

Y+,MAX  =  15.7916  N  =  150  PRINT  INT  *  2  DY+  =  1.0528E-01 


Y+  U+  (U+  -  Y+)/U+  (U+  -  Y+)/Y+  (U+  -  UC)/U+ 

0.000  0.000  0. 0000E+00  0. OOOOE+OO 

15.792  12.076  -3.0769E-01  -2.3529E-01 

ITERATION:  y ,meas( cm) ,  u,meas( cm/sec ) ,  u,tau(cm/sec) 
1.201900000000000  9.579000000000000  0.7841232469575866 

y,+(pred),  u+(pred),  UIC 

15.79159467317374  12.07598699939817  9.469062136185694 

Convergence  error( fraction)  -  -1 . 1476966678599649E-02 


MEAN  VELOCITY  PROFILE 

VAN  DRIEST  (HUFFMAN/BRADSHAW)  MODEL 

KAPPA  -  0.400000  Y+,L  -  79.0784  K,P  =  -2.0414E-02 

Y+,MAX  *  15.9025  N  -  150  PRINT  INT  «  2  DY+  -  1.0602E-01 

Y+  U+  (U+  -  Y+)/U+  (U+  -  Y+)/Y+  (U+  -  UC)/U+ 
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Typical  output  from  interactive  session  (Cont'd.) 


0.000  0.000  0 . 0000E+00  0.0000E+00 

15.902  12.112  -3 . 1293E-01  -2.3834E-01 

ITERATION:  y,meas(cm),  u,meas(cm/sec) ,  u, tau( cm/sec) 
1.201900000000000  9.579000000000000  0.7896296084287344 

y,+(pred),  u+(pred),  UIC 

15.90248824610853  12.11222110904989  9.564168411541314 

Convergence  error ( fraction)  <=  -1 . 5483441339060719E-03 


MEAN  VELOCITY  PROFILE 
VAN  DRIEST  (HUFFMAN/BRADSHAW)  MODEL 
KAPPA  =  0.400000  Y+,L  ■=  78.6299  K,P  -=  -2.0347E-02 


Y+ , MAX  = 

15.9197 

N  *=  150  PRINT  INT  -  2 

DY+  -  1 . 0613E-01 

Y+ 

U+ 

(U+  -  Y+)/U+ 

(U+  -  Y+)/Y+ 

(U+  -  UC)/U+ 

0.000 

0.000 

0.0000E+00 

0.0000E+00 

15.920 

12.119 

-3 . 1366E-01 

-2 . 3877E-01 

ITERATION:  y,meas(cm),  u,meas( cm/sec ) ,  u, tau( cm/sec ) 
1.201900000000000  9.579000000000000  0.7904848547073068 

y,+(pred),  u+(pred),  UIC 

15.91971219990579  12.11859558323355  9.579566268868983 

Convergence  error ( fraction)  =  5 .9115656016590674E-05 

MEAN  VELOCITY  PROFILE 
VAN  DRIEST  (HUFFMAN/BRADSHAW)  MODEL 
KAPPA  =  0.400000  Y+ , L  -  78.6462  K,P  -  -2.0350E-02 


Y+ , MAX  = 

15.9191 

N  •=  150  PRINT  INT  =  2 

DY+  =  1.0613E-01 

Y+ 

U+ 

(U+  -  Y+ )/U+ 

(U+  -  Y+ )/Y+ 

(U+  -  UC)/U+ 

0.000 

0.000 

0 . 0000E+00 

0.0000E+00 

15.919 

12.119 

-3 . 1358E-01 

-2 . 3872E-01 

ITERATION:  y ,meas( cm) ,  u ,meas( cm/sec ) ,  u,tau(cm/sec) 
1.201900000000000  9.579000000000000  0.7904534187114223 

y,+(pred),  u+(pred),  UIC 

15.91907910490817  12.11889022361206  9.579418208242582 

Convergence  error ( fraction)  -  4 . 3658862363722366E-05 


MEAN  VELOCITY  PROFILE 

VAN  DRIEST  (HUFFMAN/BRADSHAW)  MODEL 

KAPPA  *  0.400000  Y+,L  =  78.6924  K,P  =  -2.0357E-02 

Y+  ,MAX  -  15.9173  N  -  150  PRINT  INT  -  2  DY+  *=  1.0612E-01 


Y+  U+  (U+  -  Y+)/U+  (U+  -  Y+ )/Y+  (U+  -  UC)/U+ 
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Typical  output  from  interactive  session  (Cont'd.) 


0.000  0.000  0.0000E+00  0.0000E+00 

15.917  12.118  -3.1347E-01  -2.3866E-01 

ITERATION:  y,meas(cm),  u,meas( cm/sec ) ,  u, tau( cm/sec ) 
1.201900000000000  9.579000000000000  0.7903646186512759 

y,+(pred),  u+(pred),  UIC 

15.91729074502697  12.11849922037130  9.578033014934551 

Convergence  error ( fraction)  =  -1 . 0094843568733802E-04 


KEAN  VELOCITY  PROFILE 


VAN  DRIEST  (HUFFMAN/BRADSHAW)  MODEL 

KAPPA  =  0.400000  Y+ , L  =  78.6602  K,P  -  -2.0352E-02 

Y+,MAX  =  15.9185  N  *  150  PRINT  I NT  *  2  DY+  =  1.0612E-01 


Y+  U+  (U+  -  Y  +  )/U+  (U+  -  Y+)/Y+  (U+  -  UC)/U+ 

0.000  0.000  0 . OOOOE+OO  0 . 0000E+00 

15.919  12.119  -3.1355E-01  -2.3870E-01 

ITERATION:  y,meas(cm),  ufmeay( cm/sec ) ,  u,tau(cm/sec) 
1.201900000000000  9.579000000000000  0.7904266098460991 

y,+(pred) ,  u+(pred),  UIC 

15.91853919649905  12.11872123327992  9.578959740091384 

Convergence  error ( fraction)  *  -4.202934399’’7^3896E-06 


ICOR  =  1,  u,tau( cm/sec)  *  0.7904293028319327 

C,f (based  on  u,cl)  =  5. 3722896225529610E-03,  K,p  -  -2 . 0351766221567797E-02 

(ICOR  =  1  =  acceleration  correction  included 
=  0  -  not  included  ) 
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